SPhT-T03/187 



Renormalization of ^-derivable approximations in scalar field 

theories 

Jean-Paul Blaizot, Edmond lancu and Urko ReinossQ 
Service de Physique Theorique, CEA/DSM/SPhT, 
91191 Gif-sur-Yvette Cedex, France. 
(Dated: February 1, 2008) 

Abstract 

We discuss the renormalization of ^-derivable approximations for scalar field theories. In such 
approximations, the self-energy is obtained as the solution of a self-consistent equation which 
effectively resums infinite subsets of diagrams of perturbation theory. We show that a consistent 
renormalization can be carried out, and we provide an explicit construction of the counterterms 
needed to eliminate the sub divergences. The counterterms are calculated from the solution 
of an auxiliary gap equation which determines the leading asymptotic part of the self-energy. 
This auxiliary gap equation may be chosen as the gap equation of the massless theory at zero 
temperature. We verify explicitly that the counterterms determined at zero temperature are 
sufficient to eliminate the divergences which occur in finite temperature calculations. 
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I. INTRODUCTION 



Self-consistent approximations, such as the ^-derivable approximations, offer a way to 
deal with physical situations where the quasiparticle picture is a good starting point. They 
allow one to include large effects of the interactions into the propagators of quasiparticles, 
leaving relatively weak residual interactions among them. Such approximations were 
introduced many years ago in the context of the non-relativistic many body problem 

13,3, 

and have been reformulated later in the language of relativistic field theory 
^. They have been used recently in various field theoretical problems. In particular, 
they form the basis for a quasiparticle description of the equilibrium thermodynamics 
of the quark-gluon plasma 3, which exploits the remarkable simplifications they lead 
to the calculation of the ent.opy B (along =i„tilat hue. mod* in^p.ed by the 2- 
loop ^-derivable approximation for QCD have also been derived \^). They constitute a 
framework for variational approximations ||8| and they provide a more general alternative 
to screened perturbation theory 3|. They have been recently applied to the dynamics 
of quantum fields out of equi, 
"conserving" approximations 



Lagrangian 



3rium Q]. ^-derivable approximations are known to be 
^ , that is, they are consistent with global symmetries of the 
ll| . However since they lead to modifications of the propagators while leaving 



the vertices unaffected, they violate the Ward identities associated to local symmetries; 
consequences of such violations have been recently explored quantitatively , and shown 
to be suppressed with respect to naive estimates based on power counting. 

One of the main obstacles for implementing such non-perturbative approximations in 
quantum field theory is the difficulty of their renormalization. The central equation that 
one has to solve is a self-consistent Dyson equation, or gap equation, for the propagator. 
Such a gap equation effectively resums an infinite set of Feynman diagrams contributing 
to the self-energy, with arbitrary powers in the coupling constant. We loose the expansion 
in powers in the coupling as a tool to organize the divergences of the various diagrams. 
Indeed, at each order in an expansion of the coupling, only a subset of the diagrams 
of perturbation theory contributing to this order is effectively taken into account by 
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the gap equation. The standard proofs of renormahzabihty do not therefore immediately 
guarantee that the divergences can be systematically eliminated. Thus, with the exception 
of the simplest self-consistent mean field approximations, which can be solved easily (see 
for instance jisl. 11^). the renormalization of ^-derivable approximations has remained a 
difficult and unsolved problem. 

One of the motivations for studying ^-derivable approximations is their potential use- 
fulness in quantum field theory at finite temperature. A problem that one meets there 
is the possible occurence of temperature dependent infinities. While, on general grounds 
(see for instance one does not expect the temperature to generate any new infinities, 
it is not always straightforward to guarantee in practical calculations beyond perturbation 
theory that potentially divergent terms dependent on the temperature do actually cancel. 
This is especially the case in situations where the various divergences cannot be classified 
according to some expansion in a small parameter. (For a simple model calculation where 
one can verify this cancellation explicitly in the first non trivial order in a expansion, 
see ^sl-) It may happen then that one is led to introduce counterterms which depend on 
the temperature (see for instance Q), which obscures the physical interpretation of the 
results. 

A maior progress on this issue was achieved recently by van Hees and Knoll in a series 



of papers 
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17| is based on an expansion 



171 lisj l . The strategy that they put forward 
of the self-consistent propagator around the corresponding propagator in the vacuum. 
They were able to prove the elimination of the temperature dependent divergences, and 
have established a well defined scheme to solve the gap equation . However, sev- 

eral features of their approach remain somewhat unsatisfactory. The real time formalism 
that they use, which appears to be essential at some stages of their derivation, does not 
allow an easy comparison with more conventional field theoretical techniques. More im- 
portantly, by focusing on the temperature-dependent divergences alone, they introduce 
disymmetrical treatments of the vacuum sector and the finite temperature one. This hides 
the fact that the particular structure of sub divergences that one needs to deal with to 
guarantee the elimination of temperature dependent infinities is already present in the 
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vacuum. Thus, they do not address the issue of renormahzabihty in the vacuum sector 
itself, but take it for granted. Finally, the approach by van Hees and Knoll cannot be 
applied as it stands to massless theories, since an expansion around the vacuum (mass- 
less) propagator would be afflicted with infrared divergences. This spoils a main benefit 
of the self-consistent approximations, which is to regulate the infrared behaviour via a 
resummation of the medium effects. 

In this paper, we present a general discussion of the renormalization of ^-derivable 
approximations for scalar field theories. We do not rely on a separation of the vac- 
uum sector from the finite temperature one, and shall use the imaginary time formalism, 
making the connection with standard field theory transparent. In fact we carry out the 
renormalization at zero temperature, and proceed as in ordinary perturbation theory, by 
analyzing the sub divergences of the diagrams which are generated as one iterates the gap 
equation. The analysis and elimination of divergences is carried out according to the 



Bogoliubov-Parasiuk-Hepp (BPH) procedure [l^ by first removing sub divergences 
using previously determined counterterms and then absorbing the remaining divergences 
into corrections to these counterterms. As we shall see, vertex divergences occur, which 
are identified with those of a special Bethe-Salpeter equation for the 4-point function. The 
counterterms will be explicitly constructed, as well as a finite gap equation. At the end 
of our analyzis, it will become straightforward to show that the temperature-dependent 
contributions are finite, as expected: potential temperature dependent divergences cancel 
out as one eliminates the sub divergences. But since carried out already in the vacuum, 
our proof of the renormahzabihty of the <l>-derivable approximations is a general one, and 
applies to any of their potential applications, in particular to off-equilibrium problems. 

The outline of the paper is the following. In the next section, we recall the essential 
features of ^-derivable approximations for the 0^ scalar field theory, their diagrammatic 
content and basic definitions concerning their renormalization. In the next two sections, 
we discuss thoroughly the renormalization of simple approximations based on the 2-loop 
and 3-loop skeletons, respectively. The discussion is carried out in such a way as to 
prepare for the generalization presented in the last section. A preliminary account of this 
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work was given in [2| 



II. GENERAL FORMALISM 

As a generic example of a scalar field theory, we shall consider that described by the 
Lagrangian 

C=\{dc^f-'^ct>^-^ct>\ (III) 

and discuss only the symmetric phase in which the field expectation value vanishes, that 
is, we shall not treat the case of spontaneously broken symmetries. 



A. Skeleton expansion for the t her mo dynamical potential 

The starting point of ^-derivable approximations is the following expression for the 
thermodynamic potential f2, considered as a functional of the full propagator D 

(3n[D] = ^Tr \ogD-^ - iTrnD + , (II.2) 

where Tr denotes the trace in configuration space, and (3 = 1/T, with T the temperature. 
The trace over configuration space involves integration over imaginary time and over 
spatial coordinates. Alternatively, these can be turned into summations over Matsubara 



frequencies and integrations over spatial momenta 

q3„ 



21 



dr J d^x ^ J [dfc], (II.3) 

where V is the spatial volume, = {iun, k) and Un = rnrT, with n an even integer. We 
have introduced a condensed notation for the measure of the loop integrals (i.e., the sum 
over the Matsubara frequencies uJn and the integral over the spatial momentum k): 

n,even ^ ' 



5 



At zero temperature, /3 — oo, the loop integrals become ordinary Euclidean integrals in 
four dimensions: 

/ 1*1 -Iwr- 

and Q becomes the ground state energy. In the following we shall often use the short- 
hand notation to denote either the sum integral ()II.4|) at finite temperature, or the 
Euclidean integral pi.5|l . depending upon the context. 

The sum-integrals in equations like Eq. pi.2|l contain ultraviolet divergences and re- 
quire regularization. Unless explicitly stated otherwise, we shall assume throughout di- 
mensional regularization. 




FIG. 1: The skeleton diagrams contributing to the thermodynamical potential to order four 
loops in scalar cp'^ theory, and, below, the corresponding contributions to the self energy to order 
three loops. 

The self-energy 11 in Eq. ()II. 2|) is related to D by Dyson's equation 

D-^ = + n, (n.6) 

where Dq denotes the free propagator. 

The quantity $[-D] is the sum of the 2-particle- irreducible "skeleton" diagrams, i.e., 
diagrams that cannot be separated into two disconnected parts by cutting two lines. In 
short, skeletons are diagrams of perturbation theory in which one cannot identify self- 
energy insertions. The skeletons contributing to the thermodynamic potential up to four 
loops and, respectively, to the self-energy up to three loops, are displayed in Fig. [H In 
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Eq. ()II.2|) . ^[D] is calculated by associating full propagators D to the lines of the skeletons 
(while skeleton diagrams in perturbation theory would be calculated with free propagators 
Do). 

There are several ways to obtain Eq. pi.2|) : by integration over the coupling constant 
by a Legendre transform with respect to a bilocal source 0,0], or via a simple 
diagrammatic analysis P]. We briefly recall here a derivation based on the latter approach. 
This will give us the oportunity to introduce concepts that will be useful later on. 

Consider then the self-energy 11. From its definition in Eq. pi.6|) . it is the sum of all 

the one-line-irreducible diagrams contributing to the propagator. Clearly, the self-energy 

skeletons can be obtained from the skeletons of fl by cutting one line of the latter (see 

Fig. n for examples). Thus we have: 

mOo] 1. 



-MDo] 



(11.7) 
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where ns[-Do] denotes the sum of self-energy skeletons evaluated with free propagators. 
Now, it is easy to verify that all the diagrams of perturbation theory contributing to 11 
can be obtained from the self-energy skeletons by replacing their internal lines by full 
propagators. In other words: 



U = Us[D]=2 




(11.8) 



FIG. 2: A typical diagram contributing to the thermodynamical potential 

This simple procedure does not hold, however, for the diagrams contributing to the 
thermodynamical potential because, in the absence of external lines, the diagrams ob- 
tained by replacing the free propagators in a skeleton by full propagators are overcounted. 



Consider, as an example, the diagram in Fig. |21 One can identify on it five ^-skeletons 
(by replacing successively by propagators the parts which can be regarded as self-energy 
insertions); these skeletons belong to the set drawn in the first line of Fig. ^ The lines 
labelled a, b, and c in Fig. |21are said to belong to cycles (there are three cycles in Fig. |2I). 
The cycles have the following property: By cutting two lines of a cycle, one cuts the 
diagram into two disconnected parts. Cutting two lines belonging to two different cycles 
leaves the diagram connected. Clearly, two cycles of lines have in common at most one $- 
skeleton (identified by replacing each cycle attached to it by a propagator), but a skeleton 
may have several cycles attached to it. In other words the cycles generate a tree structure, 
whose branches are chains of skeletons attached by cycles of lines (an illustration of such 
a structure is given in Fig. E] below). This particular structure underlies the algebraic 
derivation of Eq. (III.2I) in terms of a Legendre transform with respect to a bilocal source 

For a given diagram 7, let ric be the number of cycles, Ug the number of ^-skeletons 
and rii the number of lines belonging to cycles. The following identity is easily proved by 
induction: 

nc + ng-ni = l. (II.9) 

For the diagram in Fig. |2l we have ric = 3, Us = 5, ni = 7. Let now ^(7) be the 
contribution of the diagramm 7, calculated with the rules of perturbation theory (i.e., 
with free propagators). We have the obvious identity: 

^^(7) = J2K{l)+ns{^)-nii^)]w{^). (11.10) 

7 7 

The usefulness of this relation is that the three terms in its right hand side correspond to 
the three terms in Eq. pi.2|) . as we now show 
We have : 

^d] = J2MiM7)- (11.11) 

7 

This states simply that when one evaluates the skeletons with full propagators one counts 
each diagram 7 as many times as one can identify skeletons in 7. For instance the diagram 
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of Fig. |21will occur once in the 1-loop skeleton of Fig. d three times in the 2-loop skeleton, 
and once in the 3-loop skeleton. 




Consider now the sum of the family of diagrams of Fig. El corresponding to a cycle of 
lines with an arbitrary number n [n >2) of self-energy insertions. This is given by 

^ Tr log(l + DoU) - i Tr Don = ^ ^0(7)^(7)- (11.12) 

7 

Again, the right hand side follows from the fact that the diagram 7 is counted as many 
times in the left hand side as there are cycles of lines in 7. 

Finally consider one line of a cycle, and view it as part of a full propagator closed on 
a self-energy. We have then 

-^TrUD + ^TrHDo = 5^nK7V(7), (11.13) 

7 

since 7 will be counted as many times in the left hand side as it contains lines of cycles. 

By adding together the contributions obtained in Eqs. ()11.11|) . ()11.12|) and ()11.13|) we 
recover Eq. pi.2j) up to the term term (1/2) Tr log Dq"^ which is proportional to the 
thermodynamical potential in the absence of interactions. 
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B. ^-derivable approximations 



An important property of the functional Q[D] is to be stationary under variations of 
D (at fixed Dq) around the physical propagator: 

5n[D]/6D = 0. (11.14) 

It is indeed easily verified that this equation for D is equivalent to the self-consistent 
Dyson equation (commonly referred to as the gap equation): 

D-' = D,' + U[D], n[D]=2^. (11.15) 

The physical thermodynamical potential is then obtained as the value of Q[D] for D 
solution of this equation. 

Self-consistent (or variational) approximations, i.e., approximations which preserve the 
stationarity property pi.l4jl . are obtained by selecting a class of skeletons in $[-D] and 
calculating 11 from Eq. pi.8j) . Such approximations are commonly called "^-derivable" 

3 

By selecting a skeleton in $, one obtains an approximation that effectively resums, 
through the iteration of the gap equation, an infinite subset of diagrams of perturbation 
theory. An example of such a resummation is illustrated in Fig. |3]a corresponding to 
the approximation based on the 3-loop skeleton. One recognizes the topology discussed 
above for the diagram contributing to Q: a tree structure, with skeletons joined together 
by cycles of lines. Consider now the diagrams obtained by opening one line of n(P) which 
does not belong to a cycle of lines (thereby selecting one branch in the tree alluded to 
before). Let K be the momentum fiowing through the corresponding propagator Dq. The 
resulting diagram is a contribution to a 4-point function that we shall call T{P, K) (see 
Fig. 01b). The diagram may contain cycles with more than two lines, that is, lines carrying 
self-energy insertions. Let then Fq be the set of diagrams of F which do not carry such 
insertions, i.e., each cycle of Fq contains only two lines (see Fig. |31c). The diagrams of Fq 
are two-line reducible. Let us call Aq the contribution to Fq which is two-line irreducible 
(cf. Fig. HJd): Ao(P, if) is the sum of all the diagrams contributing to Vq{P,K) which, 
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K 



(a) 





K 



(d) 



(b) 



K 



(c) 



K 



FIG. 4: (a): A typical diagram contributing to 11 when the 3-loop skeleton is chosen for <I>. (b): 
A diagram contributing to T{P, K), and which is obtained by opening the line marked by a cross 
in (a), (c): The corresponding contribution to Tq{P,K). (d): The kernel Aq whose iteration 
produces the diagram (c). 



by cutting two internal lines, cannot be split into two pieces, one containing the lines 
carrying P the other the lines carrying K. Clearly, the diagrams contributing to Tq have 
the form of chains of Aq's joined together by pairs of propagators Dq (see Fig. Elc). The 
diagrams that we have identified play the role of skeletons for F and A, respectively. As 
was the case for 11, the presence of external lines makes the identification of skeletons 
unambiguous, and the complete set of diagrams contributing to F and A can be obtained 
by replacing in Fq and Aq the free propagators -Do by the full propagator D (cf. Eq. pi.8j) 
for n). 

It it not difficult to show that A can be obtained from 11 by functional differentiation: 

52$ 



{11.16) 



6D{K) 6D{K)SD{P) 
The contributions to A, up to 2-loop, are displayed in Fig. As for the diagrams of F, 
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p 



p p 



p 



FIG. 5: The skeleton diagrams contributing to the kernel of the Bethe-Salpeter equation up two 
loops 

they are obtained by solving iteratively the following Bethe-Salpeter equation (see Fig. Ej): 
T{P,K) = A{P,K)-^ I A{P,Q)D'{Q)T{Q,K) 

= A{P,K)-^ I T{P,Q)D\Q)A{Q,K). (11.17) 

J Q 

Note that the symmetry property A(P, K) = A{K, P) entails the corresponding prop- 
erty for the full 4-point function: F(P, i^') = F(ir, P). This is easily verified by solving 
Eq. pL17|l iteratively. Note also that the Bethe-Salpeter equation pi.l7|) realizes a re- 
summation in a specific channel, which one may call the t-channel: this is the channel in 
which the kernel A is 2-particle irreducible. It follows that as soon as one restricts oneself 



to a limited class of skeletons, one looses the crossing symmetry of the general 



In 



point 



function. This is a general difficulty with ^-derivable approximations (see also 

k k k k 



k k 




P P 




FIG. 6: The Bethe-Salpeter equation in the t-channel 

Although r is a function of two momenta while D is a function of a single momentum, 
the equation for F is somewhat simpler than the gap equation for D. This is because 
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the equation for 11 is a self-consistent equation: 11 enters the calculation of the dressed 
propagator which in turns is used to calculate 11. In contrast, the 4-point function F does 
not enter the determination of the kernel A of the Bethe-Salpeter equation. 

For the construction of the mass counterterm, it will be later convenient to consider 
also the derivative of the self-energy with respect to the mass. In perturbation theory, 
the mass enters only the free propagator Do, and dDo/dm? = —Dq. Thus the diagrams 
for m^dll/dm? are those obtained by inserting m? in one of the internal propagators. A 
simple analysis then shows that 



This equation can be understood as follows: Consider first those diagrams contributing to 
m?' {dll / dm?) which contain no self-energy insertions (skeletons). It is then not hard to 
see that the contribution of these diagrams is given by the r.h.s. of Eq. pi.l8j) in which D 
is replaced by Dq, and similarly F is replaced by Fq. Eq. flII.18|) then follows by noticing 
that the complete set of diagrams contributing to m?{dTl/dm?) is obtained by replacing 
in the skeleton diagrams identified above the free propagators Dq by the full ones D. 

Eq. flII.18jl may be also obtained by the following algebraic manipulations that will be 
used later on. Using Eq. pi.l6j) and the fact that dD /dm? = — _D^(1 + dU/dm^), we get 




(11.18) 




(11.19) 



which we can rearrange as follows 



f \5iK-Q) + U{K,Q)D\Q) 

^ Q L 




[ K{K,Q)D\Q). (n.20) 



Now, the Bethe-Salpeter equation pi.l7j) can be written as 




K{P,K) 



(11.21) 



or equivalently as 




F(p,g). 



(11.22) 
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From these, one easily obtains 



K 



6{P-K)--r{P,K)D\K) 



6{K-Q) + -A{K,Q)D\Q) 



6{P-Q). (11.23) 



Using these equations, one easily transforms Eq. ()II.20|) into Eq. ()II.18|) . Renormalization 
counterterms will slightly modify the form of this equation, as we shall see shortly. 



C. Renormalization: generalities 

In ip^ theory in four dimensions, we know from power counting that, aside from the 
thermodynamical potential itself, only the 2-point and the 4-point functions can be di- 
vergent when computed in perturbation theory. Thus, one expects the sub divergences 
in any diagram of perturbation theory to be of two types, associated respectively with 
divergences of a 2-point and a 4-point function. An important class of 4-point subdiver- 
gences in a self-energy diagram will contain as external lines those of the self-energy itself 
and those obtained by cutting a line (associated with a free propagator) which does not 
belong to a cycle of lines. The analysis done previously indicates that this is nothing but 
a contribution to V. One should therefore expect a strong interplay between the vertex 
sub divergences of the gap equation and those of the Bethe-Salpeter equation. 

The divergences of the n-point functions can be absorbed in local counterterms corre- 
sponding to a redefinition of the parameters of the lagrangian. The initial theory depends 
on three parameters, the bare mass and coupling constant, and respectively, and 
an ultraviolet cut-off. The relations between the renormalized and the bare parameters 
(denoted by a subscript B) are the usual ones: 

0B = v^</', Zm\ = m^ + 5m^, Z'^Xb = ZxX = \ + 5\ 5Z = Z-l. (11.24) 

The parameters of the renormalized theory are the renormalized mass and coupling con- 
stant, m and A respectively, and the counterterms 5Z, dm?, 5\ are supposed to be func- 
tions of m, A and the cut-off. They are such that physical observables, expressed in terms 
of renormalized quantities, have a finite limit when the cut-off goes to infinity. 
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The relation between bare and renormalized propagators reads: 



Db 



ZD 



(11.25) 



and we define bare and renormalized self-energies by 



D-s\K) = K^ + ml + IiB 



(11.26) 



Thus n and Hb are related by 



n = ZUb + K^SZ + bm^. 



(11.27) 



For the skeleton diagrams, we have the identity 



^\\b,Db\=^\Zx\D\. 



(11.28) 



This follows from the fact that to each vertex of $ are attached four propagators. Thus, 
in $, one can replace Db = ZD by D provided one replaces at the same time Ab by 
Z^Xb = Z\\. This relation pi.28|) . together with the definitions above, allows us to 
express the thermodynamical potential f2 in terms of renormalized quantities 



where IS£1 is an infinite constant, independent of the temperature, whose role is to absorb 
the global divergence which remains in once the 2-point and 4-point function subdi- 
vergences have been properly eliminated. When expressed in terms of bare quantities, 
the thermodynamical potential is stationary with respect to variations of Db^ at fixed 
bare parameters (see Eq. pi.l4|) ). One can verify that a similar property holds when is 
expressed in terms of the renormalized propagator, i.e., the expression ffTr29|) is station- 
ary with respect to variations of D at fixed renormalized parameters (and therefore fixed 
values of the counterterms) . We have also 8^/5Db = (1/2)115, and therefore: 



l3{n[D]- = - [ logD-\K)-- [ {n{K)-SZK^-Sm^)D{K) + ^Zx\,D] 




(11.29) 



6<i>[ZxX,D] _ 1 

Id ~ 2 




(11.30) 
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with Ub evaluated with the propagator D. Given the relation ()II.26|) between and 11, 
one then obtains the gap equation in the form: 

U{K) = 5ZK' + Sm' + 2^^^^. (11.31) 

oD 

Consider finally the Bethe-Salpeter equation. In terms of bare quantities it reads: 

Tb{P, K) = Ab{P, K)-^ [ Ab{P, Q)DUQ)Tb{Q, K), (11.32) 

where 

A.(P,i^)=4^^^%^. (11.33) 

' ^ 6Db{P)6Db{K) ^ ' 

In terms of renormalized quantities, this is simply 

r(P, K) = A(P, K)-^ f A(P, Q)D'{Q)T{Q, K), (11.34) 

J Q 

with 

A(P.A-)^Z=A«(P.A-)^4,?™j (11.35) 

Note that A is not a finite quantity: as we shall see, it contains coupling constant coun- 
terterms which are needed to remove the divergences in the integral equation for F. 

Related to the Bethe-Salpeter equation is the equation for dU/dm^. In terms of bare 
quantities this is Eq. pi.l8j) . For the renormalized theory we have: 

,„^^ . .'4 f MK,Q)DHQ) (l + M«) V W, (11.36) 



dm"^ 2 7q ' \ dm'^ J 

where we have used the fact that 5m? /m? can be chosen to be a constant independent of 
m^. This equation will be used to determine the mass counterterm from the requirement 
that dli/dm? be finite. 

Let us now summarize the procedure that we shall follow in order to identify properly 
the sub divergences. We shall expand the self-energy in powers of the coupling constant, 
keeping at a given order only those diagrams which are generated by iterating the gap 
equation. Then we shall perform a standard analysis of the sub divergences and show 
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that they can be ehminated order by order by appropriate local counterterms. We shall 
be helped in this analysis by the consideration of the Bethe-Salpeter equation which will 
lead to the determination of the coupling constant counterterm SX, and the equation for 
dU/dm^ which will fix the mass counterterm Sm^. The last counterterm, 6Z, associated 
with field normalization, will be obtained independently by considering an equation for 
the asymptotic part of the propagator. Once the counterterms are known, one can verify 
that the gap equation is finite, and we shall present a well defined procedure to obtain such 
an equation. One can also verify that the temperature dependent infinities go away: these 
appear typically as products of finite temperature contributions, which are finite, by zero 
temperature subdiagrams, which are infinite, but are renormalized by zero temperature 
counterterms. In summary, once the sub divergences have been properly eliminated, all 
apparent temperature dependent divergences are also eliminated. 



III. THE TWO-LOOP APPROXIMATION 



In this section, we discuss the approximation obtained by keeping in $ only the two- 
loop skeleton. The presentation is more detailed than needed for the solution of the 



problem, which is well known 



141 ]. However, many features of the renormalization 
of a self-consistent approximation already appear in this simple example and can be 
made explicit with analytical calculations. This two-loop example constitutes therefore 
an excellent preparation for the generalizations to be done in the next sections. 
The contribution of the two-loop skeleton is 

The self-energy, obtained from Eq. pi.31|) . reads then: 

X + 6X 



n 



j D{P) + Sm^ (111.2) 



2 

where D^^{P) = + iv? + 11, and 11 is in this case a constant (so that 5Z = and 
Db = D). Finally the kernel of the Bethe-Salpeter equation is 

A = Ab = A + (5A, (III.3) 
17 



and the Bethe-Salpeter equation itself reads: 

T = A-^T I^D^P), (III.4) 
with r a constant. This equation can be conveniently rewritten in the form: 

The integrals above are ultraviolet divergent. We define them through dimensional reg- 
ularisation, with d = A — 2e and fi the renormalization scale. The dimension of the 4-point 
functions A and F becomes then that of the bare coupling constant in dimensions, 
i.e., It is convenient to rescale these functions so as to keep them dimensionless: 

A — > fi^'A, r — > fi^T, Xb — >fi^'{X + SX). (III.6) 

The Bethe-Salpeter equation reads then as in Eq. pil.4j) . with a factor fi'^'^ multiplying 
the integral. As for the gap equation it becomes 

n = ^/i2^ / D{P)+5m^, (III.7) 
2 Jp 

and n has the same dimension as m^. 

The foregoing analysis will involve Feynman diagrams calculated with the free propa- 
gator: Do{P) = (P^ + m'^)"^. The needed integrals are (with = Ane''^'^ fi"^): 



i3o(m^) ^ DoiP) = \ - - \n— + 1 \ ^ + b^, (III.8) 

and 



p 



Note that Bi is dimensionless and that its divergent piece Oi/e is independent of m 
(ai = 1/I67r2). 

At finite temperature, these two integrals receive additional finite contributions, re- 
spectively Bq and bi, with 

T_ r 
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and bj = —dh]^/dm?. 

Before we proceed, let us comment on a typical difficulty that we are faced with. 
Since the divergence in Eq. ()111.2|) is a local one, one could, naively, attempt to absorb 
it in the mass counterterm 5m?. However, if one does this, one obtains a counterterm 
whose divergent part depends on the self-energy 11. This is not desirable since, at finite 
temperature, 11 depends on the temperature, and that would introduce temperature de- 
pendent counterterms in the calculation. Besides, one expects on general grounds to be 
a ble to put the mass counterterm in the form 5m? = m^C(A,e), with C independent 
of m [3]. Clearly, the H-dependent mass counterterm just defined would not have this 
property, as easily seen by taking a derivative with respect to m^. In fact the infinities 
in Eq. ()111.2|) result from sub divergences which are best exhibited by iterating the gap 
equation. And to properly identify 5X and 5m^, one needs to carefully disentangle these 
various subdivergences. 

We are therfore led to analyze the formal solutions of the gap equation pil.7|) . and 
of the Bethe-Salpeter equation pil.4|l . that are obtained by the successive iterations of 
these equations. Such iterations construct the solutions as powers series in the coupling 
constant A. These series coincide with those one would obtain by applying the rules of 
ordinary perturbation theory, and restricting the Feynman diagrams to belong to a certain 
topology. The solutions after n iterations, to be denoted respectively by n*^"-* and T^^\ 
are of order A"'^^, with each iteration bringing in a new power of A, and associated new 
divergences. The renormalization consists in showing that, given the mass and coupling 
constant counterterms determined at the previous iterations, of order A" at most, the new 
divergences can be absorbed in a new contribution to these counterterms, of order A*^"^^. 
Note that whether this can be done or not is not guaranteed a priori by the standard proof 
of perturbative renormalizability, since at any order A" only a subset of terms contributing 
at this order is included. 
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A. Iterative solution of the Bethe-Salpeter equation 



We consider first the Bethe-Salpeter equation pil.4|) in which we replace the propagator 
D by the free propagator Dq. We call Fq the corresponding 4-point function, that is: 

To = A - ^ To/i^^ Dl{P) = A - ^roi3iK). (III.ll) 

As we have just mentioned, we shall construct the counterterms by building the solution 
of Eq. pil.ll|) as a formal series in powers of the coupling constant A. This is obtained 
by iteration, keeping at each iteration k only terms which contribute to Tq up to order 
\k+i_ xhus, Fq^^ the value of Fq obtained after n iterations, is of order A""*"^, and so is 
the counterterm ^A^, that needs to be adjusted at iteration n. Clearly, we can write 

n 

[A][„+i] = A + J](5Afc. (in.l2) 

k=l 

where 5Xk is the counterterm of order A*^"*"^. The notation [■■•][„], to be used throughout, 
indicates that only the terms up to order A" are to be kept within the brakets. 





+ X 



FIG. 7: The counterterm 5\i ehminates the divergence of the one-loop contribution to the 



4-point function Fq 



(1) 



In leading order, the solution Fq^'' is simply the renormalized coupling constant 



rr = A. 



(III.13) 



The first iteration gives 



F, 



(1) 



2^1 



A + ai-A"y 



A + 5Ai-^(^ + 6 



(III.14) 
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As anticipated, there is a divergent contribution given by the one-loop diagram in Fig. [7| 
This divergence is absorbed in the counterterm 6Xi. Using minimal subtraction, one gets: 

SX, = (III.15) 

Note that we have ignored in ()III.14|) the terms ~ X6Xi and ~ 6X^ which are of order 
higher than A^. Note also that T^^ is only a part of the 4-point function that would 
be obtained in perturbation theory at order A^. This is so because, as we have already 
mentioned, the Bethe-Salpeter equation performs a resummation in only one out of three 
possible channels. The contributions of other channels are recovered only when included 
in the kernel A (see next section). 




FIG. 8: The counterterm 5X2 eliminates the global divergence which remains in the two-loop 
contribution to Fq (diagram on left) after the sub divergences have been removed by the two 
middle diagrams involving 5Xi. 



Doing the second iteration and keeping only terms up to order A'^, one obtains 



^(2) 



A-A^r 



[3] 



A + 5Ai + 5X2 - A (A + 25Xi) ^ + A^ 1 



2 



4 



(111.16) 



where the notation [■ ■■][«] is that of Eq. ()111.12|) . Note that although there is a single 
counterterm to adjust, namely 6X2, one needs to keep combinations of terms which where 
ignored at the previous iteration, such as X6X1: such terms are needed to remove the 
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sub divergences that are now present, as can be seen in the diagrams displayed in Fig. |HI 
Observe that in the cancellation of subdivergences, the term involving the product of 6Xi 
and the finite part of Bi also cancels. The global divergence that remains, and that is 
absorbed in 6X2, depends therefore only on ai. A simple calculation gives 

SX. = x(^X. (III.17) 



2e 





+ 




FIG. 9: Diagrammatic interpretation of the Bethe-Salpeter equation. 



Quite generally, at iteration n one gets 



A 



AFi"-^)^' 
° 2 



fIII.181 



[n+l] 



The term Fq" is assumed to be finite. It contains all the counterterms which eliminate 
the subdivergences that do not involve the lowest vertex in Fig. El The factor A multiply- 
ing it brings the counterterms needed to remove the subdivergences involving the lower 
vertex. After all the subdivergences have been removed by the counterterms determined 
in the (n — 1) previous iterations, there remains a global divergence which is cancelled by 
the adjustment of the contribution 5A„ to A (which affects only the first A in the r.h.s. 
of Eq. ()III.18|) ). An illustration at order is provided in Fig. ^1 
One can rewrite Eq. piI.18|) as follows: 



F, 



(n) 



A 



1 + Af 



[ro] 



n+l] 



(III.19) 



[n.+l] 



The last equality emphasizes the fact that Fq is nothing but the perturbative solution 
of order A""^^, restricted to the particular resummation of the chain of bubbles in the 
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+ 



+ 




+ 



+ 



5X, 



(a) 



(b) 



(c) 



(d) 



(2) 

FIG. 10: Contributions to Tq at order A , illustrating the elimination of the various subdiver- 
gences of the diagram (a) by counterterms determined at the first two iterations (i.e., 6Xi and 
SX2)', these leave a global divergence absorbed in 8X3. Diagrams (b,c,d) represent counterterms 
for subdivergences which do not involve the lower vertex of (a); diagrams (e) and (f) eliminate 
the subdivergences involving the lower vertex of (a). 

t-channel. In order to determine the general form of the counterterms, it is in fact more 
convenient to invert this relation: 



i 0— 



(III.20) 



[n+l] 



and exploit the fact that Tq is finite. Then, it is not hard to get 



A 



To 







6Xr 



A 



Aai 



(III.21) 



where ai and bi are defined in Eq. pil.9jl . 

At this point we note that the counterterms depend only on ai, that is, they are 
independent of the mass. It follows that the Bethe-Salpeter equation pil.ll|) written 
with the full propagator D instead of the free propagator Dq, can be made finite with the 
same counterterms as those determined for Fq. This may also be seen by observing that 
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any insertion of a self-energy in the propagators of the Bethe-Salpeter equation alters only 
the finite part of the integral in Eq. pil.4|) . 

Let then F be the solution obtained with the full propagator D. From the first of 
Eqs. ()111.21|) . and a similar equation for F (obtained by replacing simply m? by m? + 11), 
one easily gets 

i = i- - i6i(m^) = i - h.im' + n). (III.22) 
Thus, F and Fq differ indeed only by a finite part: 

- - — = ^ In (III.23) 
where we have used the explicit form of 61 (m^) given in Eq. pil.9|) . 



m 




FIG. 11: Diagrammatic representation of the two terms in the right hand side of Eq. (|III.24|) . 



B. The equation for the mass counterterm 

We now turn to the equation satisfied by 811/ dm?, whose structure, as discussed in 
Sect. II, is very close to that of the Bethe-Salpeter equation. This will allow us to 
determine the counterterm 5m^, which we assume of the form 5m? = m^C(A,e). From 
Eq. ()11.19|) (or directly from Eq. ()111.2|) ). one easily gets: 
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where Bi stands here for Bi{m? + 11). This equation has a simple diagrammatic interpre- 
tation: each term represents a tower of bubbles, with or 5m? inserted in the upper 
most bubble (see Fig. ITT|). The corresponding diagrams are easily generated by expanding 
the right hand side of this equation in powers of A. A simple analysis then reveals that 
the subdivergences associated with coupling constant counterterms cancel out; that is, the 
potential subdivergences are removed by the same counterterms as introduced in solving 
the Bethe-Salpeter equation. Such subdivergences in the towers of bubbles of Fig. ITT] are 
those which do not involve the upper loop with m? or 5m? inserted. The subdivergences 
involving the upper loops are remove by adjusting the mass counterterm 5m^ at each 
iteration. The elimination of coupling constant subdivergences can be made explicit by 
using the Bethe-Salpeter equation pil.l9|) to write: 



m 



dm^ 

where we have used the identity 



-^-B^m' 



5m 



r 



^1 



1. 



(III.25) 



(III.26) 



The mass counterterms can then be easily determined by expanding the r.h.s. of 
Eq. pil.25|) in powers of A, identifying at each order the divergent part of 5mi?. In 
leading order we get 

2 Aai 



m 



2e 



At next order, we have: 

2an« 



m 



-m 



r(i) 

' B, + 5ml 



^B, ) + 5ml 



(III.27) 



(III.28) 



dm'^ 2 

where n(i) and F^^) are the values of 11 and F obtained at the first iteration of, respectively, 
the gap equation and the Bethe-Salpeter equation. Keeping only the divergent parts of 
order A^ one gets 



5ml 



2/2 



+ A 



5ml 



div 



5ml 
m^ 



Aai 
"27 



A6i 



(III.29) 
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where we have used Eq. fllll.l4|) for r"^^). Observe the cancellation of the term involving the 
leading order mass counterterm 6ml multiplied by the finite part bi, leaving the following 
formula for 5m^: 

6ml 6ml ( ■^'^i 



(III.30) 



m^ m^ \ 2e 
6\i ai 



In the second line we have written the formula which emerges naturally if one cancels 
the sub divergences in a different order (making use of the explicit expression pil.27j) of 
5mo). This writing generalizes to high orders. This can be seen by rewriting Eq. ()III.25|1 
as follows: 



m 



~m^-Bi + -6m^. (III.31) 



(9m2 2 A 

Clearly, the divergence in the right hand side is cancelled if we choose 

6m'^ A tti 



(III.32) 



m^ 2 e 

We have obtained earlier A as a power series in the coupling constant (see Eqs. ()III.12|) 
and pil.21|) ). Similarly we can write 

n 
fc=0 

with 

"itT- ,n,.34) 

As was the case for the Bethe-Salpeter equation, it does not matter for the determi- 
nation of the divergent part of the mass counterterm whether the free propagator Dq, 



or the full propagator D is used in Eq. flIII.31jl : this is a priori clear from the diagram- 
matic structure in Fig. ITD and is obvious in Eq. pil.32|l . But Eqs. (IIII.24|1 or piI.3Hl for 
dU/dm"^, including the finite parts, hold only if the full propagator is used (as is clear 
from their derivation). 
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C. Iterative solution of the gap equation 



We turn now to the gap equation ()III.7|) . and construct iteratively its solution as a 
formal series in powers of the coupling constant A. At a given order in A, 11 contains 
two types of subdivergences: those associated with the renormalization of the coupling 
constant, and those associated with the renormalization of the mass. 

We shall verify that the counterterms 5A„ determined in solving the Bethe-Salpeter 
equation eliminate the subdivergence of the first kind, while those of the second kind are 
eliminated by the mass counterterm obtained at the end of the previous subsection. 




We start by expanding the r.h.s. of the gap equation in powers of 11: 

n = - / — + 6m^ 



+ (III.35) 



2 

k=0 



The coefficients Bk are given by 



= I ... = ^ J^- (111.36) 



/p (P2 + rn2)fc+i k\ (9(m2)^' 

The expansion ()III.35|) is represented diagrammatically in FiglT^ Note that only Bq 
and Bi are ultraviolet divergent and, as we have seen in the previous subsection, in 
minimal subtraction the counterterms 6X and Sm^ depend only on their divergent parts, 
respectively and ai. Each iteration brings in one power of A so that n'^"\ the value of 
n at iteration n, is of order A"^^. This means that at iteration n > 1 one can stop the 
expansion ()III.35|) at k = n. 
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The lowest order contribution, n'^°\ is that of the 1-loop diagram in FiglT^ 

n(°) = ^B, + = ^ (I + 60) + Sml (III.37) 

There is no subdivergence in this diagram and the global divergence is absorbed in the 
mass counterterm Stuq given in Eq. pil.27|) (recall that ao = — m^ai). 



5m2 



FIG. 13: The 1-loop contribution to II^'^^ , and the corresponding mass counterterm. 




The first iteration of the gap equation, in which we keep only the terms up to order 
A^, leads to 



[2] ^ 



A H / A^ \ A 

-Bo + Sml + — { 5Ai Bi - 5ml-Bi + 5ml. (111.38) 



nd) 



2 ^ 2 V 2 7 "2 

In the second line we have grouped the terms so as to display the cancellation of diver- 
gences: one recognizes indeed in the last four terms, of order A^, the structure already 
met in Eq. pil.28j) for dTi/drn^ (recall that the divergent part of Bq is the same as that 
of m^Bi). This is illustrated in Fig. El 

A further iteration will allow us to uncover the generic features of the iteration proce- 
dure. Keeping only terms up to order A^ on obtains: 

n(^) = (n(°))^ - ^i3in(^) - ^i3,n(°) + \ [^\,, b, + M , (111.39) 

This iteration adds to 11 a 3-loop contribution, whose corresponding diagrams are dis- 
played in FigH^l The first term in Eq. pil.39|) . quadratic in 11, is finite; all the divergences 
are localized in the terms linear in 11 and to isolate them we need to express U^^^ and 11^) 
in terms of elementary integrals. Again, as was the case in the first iteration, we shall 
find that the counterterms determined at the previous two iterations eliminate all the 
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FIG. 14: The 2-loop contribution to exhibiting two sub divergences. That associated with 

couphng constant renormaUzation (the lower loop) is eliminated by diagram (b) involving 6Xi; 
that associated with mass renormalization (the upper loop in (a)) is eliminated by diagram (c) 
involving Sijiq. The remaining global divergence is absorbed in the mass counterterm (5mf . 




FIG. 15: The 3-loop contribution to n^^-*, with its coupling constant counterterms (mass coun- 
terterms are not drawn). 

subdivergences but one, namely that associated with diagram (a). This subdivergence is 
recognized as the global divergence of the order contribution to the 4-point function 
and it is absorbed in 6X2- The remaining divergence is a global one and is eliminated by 
new mass counterterm 61712. 

It is useful to observe that the new divergences which enter at each new iteration n, 
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and which require the new counterterms 5A„ and Sm"^, are entirely contained in the tower 
of bubbles such as diagram (a) in Fig. . The other structures which develop as we 
proceed through the iterations, such as diagram (e) in Fig. E[ are all made finite by 
previously determined counterterms (for diagram (e) this is due to the mass counterterm 
that is not drawn here). Then, it is not hard to see that all the divergences in such a 
tower cancel in exactly the same way as they do in the equation for dU/dm'^. To make 
this explicit, note that we can rewrite Eq. pil.39|) as follows 



nP' = ^B2(n(»>)' + 



2/2 



J [3] 



2 



-1 



(III.40) 



[3] 



The last two terms represent the contribution to 11*^^) obtained by iterating only the term 
linear in 11, thus generating the tower of bubbles such as that in Fig. llBr a). together with 
the counterterms which make it finite. That the contribution of these last two terms of 
Eq. pil.4njl is finite is made clear by comparing them to the right hand side of pil.24j) . 
observing that factor in front of Bq in Eq. pil.40|) is r/2 (see also Eq. pil.25|l ). and that 
the divergent parts of Bq and —m'^Bi are the same. 

One can summarize the iterative procedure by the following formula giving the solution 
n*^") at iteration n in terms of the (finite) solution n*^"^^) obtained at iteration (n — 1): 



n(") 



A 



k=l 



B 



[n+l] 



+ [A][n+l]Y+M[.+l], (111.41) 



where we have separated the terms which contain n^""^^ from those which contain only the 
new counterterms that need to be introduced at each iteration, and which are independent 
of n. It will prove convenient to introduce a special notation for these two contributions. 
We set: 

^A " 



n 



(n) 



Y^i-i^B, (n("-i))' 



k=l 



J [n+l] 



n 



(n) 



Bo 

[n+l] 2 



n+l] 



(III.42) 



(III.43) 



so that n(") =llW 
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From the previous analysis, we expect the new sub divergences to occur in the term 
with k = 1. However, we have to make sure that sub divergences introduced by the factor 
A in front of otherwise finite integrals do cancel. In order to do so, we use Eq. fjlll.l8|) to 
express A in terms of Fq within Eq. ()111.42|) . and obtain: 



n 



in) 



k=l 



2 

k Fn 



[n+1] 



k=2 



(111.44) 



[n+1] 



The first term is now finite. The divergences have been isolated in the last two terms 
where we can identify nQ"""^-* — n^""^) = — fl^"'^"'^-*. One then gets for n*^"-) = Hq"^ + 112"'': 



k=2 



-I [n+1] 



2 2 



(III.45) 



[n+1] 



Using the definition (llll.43jl of Hg, one can combine the last two terms as follows 

Bo ' 



A - ^B.A 



+ 



J [n+l] 



J [n+1] 



(in.46) 



But the divergent terms in this expression are those of the right hand side of Eq. pil.25|) 
(see the discussion after Eq. (IIII.40|) ). The two terms in Eq. f)III.46|) therefore add up to a 
finite quantity, [Fo]j„^i] bo/2. We are then left with a finite equation for determining the 
finite part of 11*^"): 



k=2 



[n+1] 



+ [ro][n+i] 



(III.47) 



This completes the proof that II*^'^-' is made finite with the counterterms determined in 
the previous subsection. 



D. Algebraic solution and scheme dependence 



We now show that is is possible to obtain the previous solutions of both the Bethe- 
Salpeter equation and the gap equation by doing simple algebraic manipulations on the 
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equations themselves, without going exphcitly through the iterations. We shall call M 
the physical mass, given by = + 11. Both M and F are physical quantities that 
are calculated as a function of the renormalized mass and coupling constant, m and A, 
using the self-consistent scheme. 

First we note that the renormalized coupling constant A is related to the bare one 
\b = /U^'A by 

1 1 ai 1 1 , , 

This follows immediately form Eqs. pil.5|) and ()III.22|) . It can be also obtained directly 
from the Bethe Salpeter equation pil.ll|) . observing that this equation is finite when 
expressed in terms of the renormalized coupling constant defined in Eq. pil.48|) . Since 
A = A + 5A, Eq. ()lll.48jl is equivalent to 

<5A = ^ = (III.49) 

which corresponds indeed to the series of counterterms that we have obtained via the 
diagrammatic analysis (cf. Eq. pil.21|) . 

Using the fact that Xb = /U^'^A is independent of /x, one easily obtains that A(/i) is 
running with n according to (we take the limit e ^ 0) 

Note that this relation is "exact" in the sense that it represents faithfully the relation 
between Xb and A, to all orders in A, within the class of diagrams which are taken 
into account in the present approximation. In this respect, it is instructive to verify 
how this relation holds as one iterates the Bethe-Salpeter equation in Sect. IIII AI By 
demanding that the solution Fq"'' be independent of jj, at iteration n (at fixed m^), the 
corrections to Eq. pil.50|) are pushed to order A"^^; that is, potential corrections of 
order A""*"^ cancel out, in very much the same way as sub divergences do. Now, while 
exact in the sense specified above, the P function in Eq. pil.50|) represents only one 
third of the perturbative one-loop /5-function. This is so because, as already mentioned, 
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the resummation performed by the one-loop gap equation involves the one-loop vertex 
correction in one channel only, namely the t-channel. The remaining two thirds of the 
one-loop /3-function will be recovered after taking the 3-loop skeleton diagram into account 
(see next section). 

Written as a function of the renormalized coupling constant, the renormalized 4-point 
function takes the form 



A 



1_ A i^i^' 



(III.51) 



Thus r is nothing but the running coupling constant at the scale ft = M. One can verify 
that r is independent of if M is, which we shall verify soon. Note that A(/i) diverges 
when ln(/i^/M^) = 327r^/r. The appearance of this "Landau pole" at large scales is a 
well known feature of 0^ field theory in four dimensions. Another manifestation of the 
same difficulty is visible on Eq. pil.49|l : when e— >0, 6\ ^ —A or \b — >0, reflecting the 
so-called "triviality" of the theory. 

Knowing F, one could easily get a finite equation for dU/dm'^. This is obtained from 
Eqs. firirni) and (ITTL32ll : 

^ = -irMM^). (111.52) 

However, it is not easy to exploit this equation to obtain a finite equation for 11 itself, 
by using methods that we could export to the general case. We shall therefore base our 
construction of a finite gap equation on the iterative solution obtained in the previous 
subsection. 

There are two observations which will guide this construction. The first observation 
is that, in the iterative solution, the cancellation of ultraviolet divergences only takes 
place up to some order dictated by the number of iterations that have been completed: 
at iteration n new divergences appear, of order A""^^, and some are ignored because they 
are of order higher than A""*"^ (these high order contributions reenter the calculation as 
countertems needed to cancel sub divergences in the following iterations). Thus, in an 
explicit algebraic solution, one should expect the cancellation of divergences to occur 
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only for the solution of the gap equation, which effectively sums up all the iterations. The 
second observation concerns the particular role of the term linear in 11 in the expansion 
of the right hand side of the gap equation in powers of 11 (see Eq. ()111.35|) ). As we have 
seen, the new divergences in the gap equation at each new iteration come from diagrams 
which are obtained by iterating just this term (e.g. the towers of bubbles in Fig. 
Terms with several insertions of 11 on the same propagator are finite (assuming that such 
insertions have been made finite in previous iterations). With these observations in mind, 
we now proceed to a construction of the algebraic solution which will follow closely what 
has been done in the previous subsection. 

We start by expanding the self-consistent propagator D around the free propagator Dq 
and set (see Eq. (1111. 35 jl ): 

D{P) =Do + 6D = Do{P) + Do{P){-U)Do{P) + Dr{P) (III.53) 

where we have singled out the term linear in 11, and Dr{P) behaves as IP/P^ at large 
P, so that the integral jpDr{P) is finite. Then we generalize the definitions pil.42|) and 
fllll.4311 and set: 

flo = ^ j^6D{P), (III.54) 

and 

U2 = ^ j^DoiP) + Sm\ (III.55) 

In these notations, the gap equation reads 

n = Ho + fl2, (111.56) 

where only IIq depends on the solution 11 (112 is entirely determined by Dq, as is obvious 
from Eq. pil.55|) and the fact that the counterterms 5A (in A = A + 6X) and Sm^ are 
calculated entirely from Dq). 

Next, one uses the Bethe-Salpeter equation to eliminate A in favor of Fq in the defining 
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equation for Ilo (Eq. ()III.54|) ) and get 

A 



5D{P) 



(III.57) 



At this point we observe that when 11 is solution of the gap equation pil.56|) . Ho — 
n = —112, so that the divergence of Ho (isolated in the second term of Eq. ()111.57|) ) is 
independent of 11: this is the main achievement of the elimination of vertex sub divergences. 
The gap equation then reads: 



(111.58) 



The first term, which involves 11 (through Dr) is finite. In the second term, we note that 
the factor multiplying 112 is Eq/A: this follows from the Bethe-Salpeter with propagators 
Do. Now 



En^ 



Eo M 



-Bo + bm' 



A ^ A \2 "^ ' ' " J 2 

where bo is the finite part of Bq, and we have set (see Eq. pil.32jl ) 

^ 2 
dm 



(III.59) 



(III.60) 



The gap equation takes then the form: 
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^0 /" / T-,^ , Eo 2 

327r2 



2 jp 
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(111.61] 



which is now a finite equation. It can also be written in terms of the renormalized coupling 
constant: 



n + n / + n 

-r = TT- In 7^ 1 

A 327r2 V P'^ 



-&o(M2 



(111.62) 
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or in terms of the physical 4-point function F: 

This last form, in which all explicit /x-dependence has disappeared, makes it obvious 
that M is independent of /i, because m^/A is independent of /i. Indeed, w^q/Xb = 
{m? + 5m^)/(yU^^A) has a finite limit when e — >■ 0, equal to m?/\. Since m?B/^B is 
manifestly independent of /i, so is therefore m^/A, which further implies that m runs with 
H in the same way as A does (i.e., according to Eq. pil.50|) ). 

As made obvious by the previous discussion, the gap equation can be written either in 
terms of Fq or in terms of F. This reflects a property that we have repeatedly emphasized, 
namely that the divergent parts of the counterterms do not depend on whether D or Dq 
is used in their calculations. We can exploit this freedom in rearranging the expansion of 
Eq. ()111.53jl in the following way: 

D{P) = Do + 6D = Do{P) + D{P){-U)D{P) + D',{P) (III.64) 

where D'^{P) behaves like 1/P^, so that the integral jpD'^{P) is finite. Then the full 
propagator D will naturally appear in the Bethe-Salpeter equation for the vertex coun- 
terterms of the gap equation. We can verify that all the steps leading from Eq. ()III.53|) to 
the final gap equation pil.63|) still hold, with the appropriate replacements of Fq by F, Dq 
by D, and by D[ (note that 112 is not affected). The gap equation pil.63|l remains of 
course unaltered. The advantage of this new expansion becomes apparent when studying 
the massless limit to which we proceed in the next subsection. 

Finally, we note that we could have used a different renormalisation scheme such as 
the mass-shell renormalization, in which we assume that m? is the physical mass, i.e., 

= M^. In this scheme 11 = 0, and the mass counterterm is given by 

The coupling constant counterterm can still be chosen as in eq. pil.48jl . and we have 
F = Fo = A(/i = m). This mass-shell scheme, which makes the gap equation trivial at 
zero temperature, simplifies it at finite temperature, as we shall see in the next subsection. 
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E. Finite temperature and massless limit 



The finite-temperature version of tlie gap equation pil.2|) reads: 

n = D{P) + ^1"^+ (III.66) 

wfiere n{ep) = l/(e^^'' — 1), Ep = a/pM-~M^ and denotes a 3-dimensional momentum 
integral. Tlie first integral jpD{P) is to be understood as before as a (i-dimensional 
Euclidean integral, but D now depends on the temperature. 

One could naively fear that temperature dependent divergences could appear in the 
calculation of 11, as for instance in the second term of Eq. pil.66j) . due to the presence of 
the factor A. This however does not happen if one carefully eliminates all the subdiver- 
gences which emerge when solving the gap equation iteratively. To understand this, it is 
enough to consider the first two iterations. 

The zeroth iteration is given by Eq. pil.37j) . with Bq replaced by bo + Bq . Clearly 
the countertem pil.27|l makes 11'^°^ finite. The first iteration is given by Eq. pil.38j) . 
with the finite parts of the integrals including their finite temperature contributions. 
Potential difficulties could arise for instance from diagram (a) in Fig. El in which the 
temperature dependent piece of the upper loop multiplies the divergent lower loop. But 
it is easy to verify that this does not happen: indeed the divergence of the lower loop is 
cancelled by the coupling constant counterterm of diagram (b). Similarly, the contribution 
proportional to the mass counterterm 6mQ multiplying the finite temperature part of the 
loop integral cancels by the same mechanism as at zero temperature. Thus, for the same 
reason that once sub divergences are eliminated, there remains no finite part coming from 
some subintegrations multiplying divergent part coming from other subintegrations, there 
does not remain either finite temperature part multiplying divergent integrals or infinite 
counterterms. This analysis can be generalized, by going through the same steps as in 
subsection IIII C\ in order to verify that, iteration by iteration, all potentially temperature 
dependent infinities actually cancel. 

To obtain this result in a more explicit fashion, we observe that the term with a 
statistical factor in Eq. pil.66|) can be viewed as a new contribution to D^, which is then 
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easy to incorporate in the renormalized gap equation ()III.62|) : 

n _ m2 + n / m2 + n \ l r n{ep) 



(III.67) 



In the mass-shell renormalization scheme, we get instead 




(111.68) 



p 



which implies that 11 ^ as T ^ 0. 

Consider now the massless limit m —>■ 0. In this limit, many of the intermediate 
manipulations based on the expansion ()III.53|) of the full propagator D around the free 
propagator Dq become meaningless because of infrared problems. For instance, the Bethe- 



therefore so is also the 4-point function Tq. More generally, all the integrals Bk with 
k > 1 become infrared divergent when m = 0. In particular, one cannot take the limit 
m ^ on the gap equation written in the form pil.68|) (or pil.61|) , since both Tq and the 
integral involving Dr develop infrared divergences in this limit. Still, the gap equation has 
a well defined massless limit, as can be seen by taking this limit at the level of Eq. ()III.67|) . 
This gives: 



which shows that, at finite temperature, the finite self-energy 11 provides an infrared cut- 
off, allowing the limit m — to be taken. This limit can also be written in terms of the 
full 4-point function F: 



Here, F is given by Eq. ()III.51|) with = 11. Thus F itself depends upon the temperature. 
From either Eq. ()lll.fi9jl or Eq. (I111.7()j) one can check that, at zero temperature, 11 = in 
this massless case. As we shall see in Sect. IV, the general approach to renormalization in 
the massless case is based on the modified expansion in Eq. pil.64|) . which takes advantage 
of the presence of the finite-temperature self-energy 11 as a natural infrared cut-off. 



Salpeter equation with propagator Dq, Eq. pil.llj) . is ill-defined when m 



0, and 




(III.69) 




(III.70) 
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IV. THE THREE-LOOP SKELETON 



Many of the developments of the previous section can be carried over to the discussion 
of the three-loop skeleton. In particular, the determination of the coupling constant 
and mass counterterms will go through the solution of the Bethe-Salpeter equation and 
the equation for 811/ dm? ^ with little structural changes with respect to the two- loop 
examples. The essential new aspect of the three-loop example is the presence of the 
field renormalization constant 5Z, and the fact that the self-energy contributes to the 
asymptotic behaviour of the propagator. 

The three loop skeleton for $ is displayed in Fig. ^ and its contribution reads: 

^^~^LIq ^(^)^(^)^(^)^(^ + + ^^^0 ' ■ ^^^■^^ 

We have added here in $ a counterterm which has the structure of the two-loop skeleton; 
this term is needed to cancel sub divergences, as we shall see. Such a term would be 
automatically present if we had included in $ both the three-loop and the two-loop 
skeletons. The self-energy follows then from Eq. pi.31|) : 

UiK) = [ [ DiP)DiQ)D{K + P + Q)+ ^-^ [ D{P) + 5m^ + K^5Z, (IV.2) 

G JpJq 2 Jp 

and the kernel of the Bethe-Salpeter equation from Eqs. pi.35|) : 

A(P, K) = -\^ [ D{Q)D{K + P + Q) + 5\. (IV.3) 

jQ 

It is important to observe that the vertices of the three-loop contribution in $ do not 
carry coupling constant counterterms. We shall see indeed that, as as we proceed to 
solving the gap equation through iterations, we do not generate structures which would 
renormalize these vertices. Such structures would appear only at the 4-loop order in the 
skeleton expansion. This is a particular illustration of a general feature of ^-derivable 
approximations that will be discussed in the next section. 

As we indicated in Sect. Ull and verified explicitly in the two loop example of Sect. IIIIl 
the sub divergences associated with vertex and mass renormalizations can be determined 



39 



from the analysis of the Bethe-Salpeter equation and of the equation for the derivative of 
the self-energy with respect to m^. We shall then go through such an analysis. But first 
we discuss in some details the first orders in perturbation theory obtained from the first 
two iterations of the gap equation. 



A. Perturbation theory 

The zeroth order iteration Il^^\K) is obtained by replacing D hy Dq = (P 



2 I ^2 



in the r.h.s. of eq. pV.2|) . and coincides with the standard 2-loop contribution of pertur- 
bation theory. In dimensional regularization, the divergent part of the diagram is of the 



form 



(IV.4) 



;i67r2)2 [24e 

Hidden in these various terms there are three distinct divergences. One is a subdivergence 
which is removed by a coupling constant counterterm. It will be dealt with shortly. 
The other two divergences are global divergences corresponding respectively to a mass 
divergence and a field normalisation divergence. These will be removed by adjusting 
respectively Sm^ and 6Z. The contribution of 6Z is easy to identify since it is associated 
with the unique divergent term proportional to K^. Using minimal subtraction we get 
the following lowest-order contribution 6Zq to 6Z : 

The adjustment of Sm"^ can be made only after the proper elimination of the vertex 
subdivergence. The latter singles out one line in the sunset diagram, namely, that line 
along which the momentum is kept finite, while the other loop momentum is sent to 
infinity. In fact, the corresponding subgraph can be identified with Aq (see Fig. Ilfij) . the 
leading order contribution to A, to within the counterterm 5A; that is: 

^fPT<^^ AX M ['a 1 |m2 + 3:(l-a;)(P + ir)2| 

Kq{P,K) - 5X = --—^ + i dxln — . (IV.6) 
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FIG. 16: The dotted line box encloses the subgraph whose divergence is that of the 4-point 
function Aq. 



Let us call 6Xq the value of 6X that renders Aq finite. Minimal subtraction gives: 



Note that this contribution, when added to that of Eq. pil.48|) . reconstructs the full order 
counterterm (3A^/32'7r^e). The counterterm pV.7j) generates an order contribution 
to n of the form: 



which eliminates the sub divergence in 11, Eq. pV.4|) . and the logarithmic term there (via 
((5Ao/2)6o involving the finite part of Bq in Eq. pV.8|) ). There remains then only a global, 
mass divergence, which is removed by appropriately tuning the counterterm 5m^. Minimal 
subtraction gives: 



Note that neither 5Ao, SZq nor drn^/m^ depend on m: these counterterms depend 
only on the asymptotic form of the propagator. Note also that 5Zq could be obtained by 
analyzing the divergence which remains after differentiating n^^-' with respect to K^. Since 
the coupling constant counterterm is independent of momentum, it does not contribute 
to the latter derivative, and the same applies to the mass counterterm. This explains why 
the determination of 5Zq can be made independently of the elimination of the vertex and 
mass subdivergences. This property persists in higher orders, as we shall see in Sect. II V CI 




(IV.8) 




(IV.9) 
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(a) 



(b) 




(d) 




6n^ + k 6Zq 




(c) 



6m^ + k^6ZQ 




FIG. 17: The diagrams contributing to the solution to the gap equation to order A'*. 

Diagram (a) is the original diagram with one insertion. All other diagrams involve counterterms 
which remove the various subdivergences. 

For later purposes, we note that the asymptotic behaviour of Aq (after renormalization) 
is of the form An(P, K) ~ ln(P + i^)^/ /i^ at large momenta, in agreement with Weinberg's 



theorem 



2J] (that Ao(P, K) is a function oi P + K can be seen directly on Eqs. pV.6|l 



or ()1V.3|) ). As for Il^^\K), its dominant contribution at large K goes as K^lnK^ (this 
may be obtained by calculating the finite part of U^^^K) for m = 0). 

Consider now the first iteration in which we keep only the contributions of order up 
to A^, and call n*^^-* the corresponding value of 11. The associated diagrams are shown 
in Fig. El Not represented there are the diagrams involving the new counterterms that 
need to be adjusted at this iteration, namely, 6Xi, 5m\^ and 5Zi. 

Let us focus on the new vertex subdivergences. The corresponding diagrams are dis- 
played in Fig. ^1 The subgraph within the box in Fig. E|a is logarithmically divergent 
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because the inserted self-energy grows like In K'^ at large as mentioned earlier; 
thus, this insertion does not change the degree of divergence of the diagram, the factor 
i^^lni^^ cancelling the l/K"^ coming from the extra propagator. The subdivergence in 
Fig. ^|a (together with the corresponding counterterms from Figs. \n\h and^lc) can 
be recognized as the contribution of order to the (global) divergence of the integral in 
Eq. ()IV.3|) for the kernel A of the Bethe-Salpeter equation. In other terms, the subgraph 
singled out in Fig. ^|a represents the first correction to the zeroth order kernel Aq in 
Eq. pV.6|) . and its divergence can be absorbed in a correction ^Aq^'' of order A^ to 5Ao. 
This type of correction is a new feature which did not appear in the two-loop example: it 
is related to the fact that the insertion of a self-energy affects the asymptotic behaviour of 
the propagator. However, the other 4-point divergence, corresponding to the subgraph in 




(a) (b) (c) 



FIG. 18: Vertex sub divergences in the diagrams contributing to 11^^^. The subdivergence repre- 
sented by the dashed- line box in figure (b) is that of the 4-point function depicted in figure (c), 
and is removed by a counterterm determined by the Bethe-Salpeter equation. 

the box in Fig. ^jb is analogous to those encountered in the two-loop example. It is the 
divergence of the 4-point function shown in Fig. ^Ic, obtained as the first iteration of the 
Bethe-Salpeter equation with the leading order kernel Aq. The corresponding divergence 
is absorbed in the counterterm 5\i (after sub divergences are eliminated by counterterms 
such as those in Figs. HZIb, d and e). 

To summarize, at the first iteration, and order A^, the vertex counterterm receives two 
contributions: one, which we call 5Aq^\ removes the order-A"^ divergence from the kernel 
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A, the other, which we call SXi, renormalizes the first iteration of the Bethe-Salpeter 
equation. Once the vertex sub divergences are removed, the mass counterterm 6m1 can 
be chosen so as to remove the remaining mass divergences from 11'^^) . Finally, SZi can be 
obtained by analyzing the divergence which remains after differentiating H^^^ with respect 
to K^. The same argument that we used at leading order shows that this can be computed 
independently of the vertex and mass renormalizations, that is, 6Zi can be obtained by 
extracting the divergence proportional to K"^ in the sum of the two diagrams (a) and (c) 
of Fig. El (with only K'^6Zq in diagam (c)). 

In order to obtain a more systematic determination of the coupling constant and mass 
counterterms, we turn now to the analysis of the Bethe-Salpeter equation and the equation 
for dU/dm'^. We shall return to 6Z in Sec. II V CI 

B. The Bethe-Salpeter equation 

As we did in the previous section, we shall consider the solution obtained by iterations, 
using first in the integral of Eq. ()II.17|) the free propagator Dq. We shall, as before, call Fq 
the corresponding solution, Aq the associated kernel, and construct Fq as a formal series 
by solving the Bethe-Salpeter equation by iterations. We have: 

Fo(P,ir) = Ao(P,K)-^ f A,iP,Q)DliQ)ToiQ,K) 

J Q 

= Ao{P,K)-^ I To{P,Q)Dl{Q)AoiQ,K). (IV.IO) 

In leading order, and before renormalization, Fq"'' is equal to Aq, which is divergent. 
This divergence has been dealt with in Sect. IIV Al where we have determined the coun- 
terterm 6X0 needed to absorb it {6X0 is given by Eq. (|IV7j) ). We shall call Aq the finite 
quantity 

Ao{P,K) = SXo-X^ [ Do{Q)Do{K + P + Q). (IV.ll) 
Jq 

Thus, after renormalization, Fq°'' = Aq. The quantities Aq and Aq differ by the counterterm 
difference 6X — 6X0 which will be needed to remove further divergences appearing as we 
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proceed with the iterative solution of Eq. ()IV.10|1 . 




FIG. 19: The contributions to Tq obtained from the zeroth and first iteration of the Bethe- 
Salpeter equation. The horizontal double bar represents Aq. The counterterm 5Xi absorbs the 
divergence of the loop integral in the last diagram for Fq^^ (see Eq. H1V.12() 1. 

The first iteration reads 

T^o\P,K) = Ao{P,K) + S\^ - ^ f K,{P,Q)Dl{Q)l^{Q,K). (IV.12) 

J Q 

The right hand side contains no subdivergences: these are ehminated by the counterterm 
which makes Aq finite. What remains in the integral is then a global divergence which 
can be absorbed in 6Xi. That the divergence of the integral in Eq. ()IV.12|) is a local 
one can be easily verified by observing that the derivative of this integral with respect to 
leaves a finite integral. This can be checked by using the explicit form of Ao{P,K) 
given in Eq. ()1V.11|) . This implies that the counterterm cannot depend on P. For the 
same reason, it cannot depend on K. Hence it is a constant, which can be chosen by 
minimal subtraction, as we did earlier for 6X0, or by imposing a renomalization condition 
on Fo(P, -ft') at some values of the momenta. We shall implicitly continue using minimal 
subtraction. A diagrammatic interpretation of the first two iterations is given in Fig. ^| 
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It is instructive to perform another iteration and get: 

Ao{P,Q) + 6X,]dI{Q) Ao{Q,K) + 6X, 



C(P,ir) = Ao{P,K) + 6X 
1 
2 
1 



Q 



+ 




Ao(P, Q)D',{Q)Ao{Q, R)D',{R)Ao{R, K), 



(IV.13) 



Q JR 



where we have dropped terms of order or higher. The contributions of order are 
depicted in Fig. 1201 One can then again proceed to an analysis of the subdivergences. Note 
that there are no subdivergence involving only one line of a Aq (since opening one line of a 
Aq would lead to a convergent 6-point function); in other words the divergent subdiagrams 
contain complete Aq as subdiagrams [l^. Again, one finds that all the subdivergences are 
eliminated by the counterterms determined at the previous iteration, namely ^Aq and 5Ai. 
The remaining global divergence is absorbed in the new local counterterm 8X2- That this 
divergence is a local one can be verified by the same argument as given after Eq. pV.12|) . 
namely by differentiating with respect to the external momenta and verifying that all the 
resulting contributions are finite. 




fi (2) ■ 

FIG. 20: The contributions of order A to Fq obtained from the second iteration of the Bethe- 
Salpeter equation, illustrating the eUmination of the subdivergences leading to the determination 
of the global counterterm 8X2- The kernel Aq (which includes its counterterm 5Xq) is represented 
by a horizontal bar, as in Fig. 1191 The black dots represents the counterterm (5Ai, except for 
the first one on the left which is the sum of 5Xi and 5X2- 



It is clear that the successive iterations of the Bethe-Salpeter equation will reproduce 
a structure analogous to that obtained in Sect. II 1 1 CI for the two-loop skeleton, with the 
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contribution Fq"^ at iteration n being of order A^""*"^. In particular, the contributions 
involving the product of a counterterm 6Xk by a finite integral cancel out (such terms 
appear for instance in the second line of Eq. pV.13|) and compensate contributions of 
the first line where one of the integration variables is kept fixed while the other one is 



taken to infinity). This is the kind of cancellations expected in the BPH analysis [la], and 
guarantees, in particular, that when the calculation is done at finite temperature, there 
does not remain, at any iteration, terms involving the product of an infinite counterterm 
by a finite temperature contribution. 

The only significant difference with the 2-loop example discussed in the previous sec- 
tion, is that the divergences do not simply factorize because the kernel is momentum 
dependent, and no simple expression can be given for the counterterms. Because of this, 
one may find useful to have finite equations for F in which those counterterms do not 
explicitly enter. Consider for instance Fo(0, P). This can be obtained from the following 
finite equation: 



Fo(0, P) - Fo(0, 0) = Ao(0, P) - Ao(0, 0) 

1 

~~ 2 



Fo(0, Q)DliQ){A,iQ, P) -Ao(Q, 0) }. (IV.14) 



It is easily verified that Ao(Q, P) - Ao{Q, 0) ~ P/Q for ^ p2_ ^^^^^ Yq{0, Q) grows at 
most logarithmically (this follows from Weinberg's theorem Q,I3), the integral over Q 
in Eq. ()IV.14|) is indeed finite. Thus, Eq. pV.14|) is finite, and can be solved to determine 
Fo(P, 0) in terms of Fo(0, 0), which is supposed to be finite and may be chosen to satisfy a 
renormalization condition. Once Fo(0,P) is known, it may be used to determine the full 
momentum dependence of Vq{K,P). The difference Vq{K,P) — Fo(0,P) is indeed given 
by the finite equation 

Fo(ir, P) - Fo(0, P) = Ao(ir, P) - Ao(0, P) 

^ ^{Ao(K,g)-Ao(o,g)}D2(g)Fo(g,p). (iv.is) 
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By combining Eqs. ()IV.14|) and ()IV.15|1 . one gets also: 

To{K, P) - ro(0, 0) = Ao{K, P) - Ao(0, 0) 

/" {Ao{K,Q) -Ao{0, Q)}Dl{Q)To{Q,P) 

-\l ro(0, Q)Dl{Q) {Ao(g, P) - Ao(g, 0)} , (IV.16) 

n 

an equation similar to that first obtained by van Hees and Knoll |17| . 

At this point, we note that the previous analysis departs somewhat from strict pertur- 
bation theory: we have indeed left out at each iteration perturbative contributions which 
correspond to self-energy insertions in the propagators of the Bethe-Salpeter equation, 
or in the lines of Aq. Let us then consider the effect of such self-energy insertions. We 
assume that the added self-energy contains all the counterterms needed to make it finite. 
The question arises then of whether this insertion brings in new types of sub divergences. 
The answer is clearly negative. Indeed, a new subdivergence could be obtained only by 
keeping finite the momentum in one of the internal line of the added self-energy. This is 
equivalent to opening that line, leading to a convergent 6-point function. It follows that 
the structure of sub divergences exhibited in the previous analysis is not modified by the 
replacement of Dq by another propagator D. However, contrary to what happens in the 
case of the 2-loop skeleton where a similar insertion of a self-energy on the intermediate 
propagator of the Bethe-Salpeter equation only affects the finite part of F, here the di- 
vergent part is also modified. This is because the dominant asymptotic behaviour of D is 
altered when the two-loop self-energy is included in the propagator. This implies that the 
numerical values of the counterterms obtained at a given iteration of the Bethe-Salpeter 
equation, such as 5Xq in Eq. pV.ll|) . or 5\i in Eq. pV.12|) . will receive corrections to 
all orders in the coupling constant coming from successive insertions of the type we have 
just discussed. We have already met an example of such a correction in the previous 
subsection, namely the correction to 5\q (cf. Fig. ITHla). 

It is then useful here to return to the gap equation in order to see diagrammatically how 
these various corrections show up as we proceed through the iterations of this equation. 
For illustration, some typical diagrams contributing to 11*^^-' at order are displayed in 
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Fig. ED The vertex sub divergences are the global divergences of the 4-point subgraphs 
enclosed within the dashed line boxes. Those in Figs. |^a-c correspond to a renormaliza- 
tion of the kernel A of the Bethe-Salpeter equation: they can be absorbed in a correction 
6X^Q^ of order to 6Xo- The two boxes in Figs. |^d and e correspond to propagator 
corrections on the first iteration of the Bethe-Salpeter equation: they contribute to the 
correction 6X^i^ of 6Xi (compare to Fig. ^|b and c). Finally, the subdivergence in Fig. 
I^ f corresponds to the second iteration of the Bethe-Salpeter equation: the corresponding 
divergence is aborbed into 6X2 (after proper elimination of the subdivergences). 






(a) 



(b) 



(c) 






(d) 



(e) 



(f) 



FIG. 21: Diagrams illustrating some of the new subdivergences which occur at order A^. The 
divergent subgraphs are enclosed in the dashed line boxes. 



We turn now to the mass counterterm Srn^. As we have seen, this can be deduced 
from the equation obtained by differentiating the gap equation with respect to m^. We 
can rewrite Eq. pi.36|) as follows: 



Q 



6{K-Q) + -A{K,Q)D\Q) 



dU{Q) 
dim? 
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A{K,Q)D\Q) + —-. (IV.17) 



Q 



This equation has a simple diagrammatic interpretation similar to that given in Fig. El 
with now A replacing the coupling constant of Fig. ^2 The diagrams contributing to 
811/ dm? that are generated by this equation have vertex and mass subdivergences, anal- 
ogous to those of the self-energy itself. It is however easy to verify that the vertex 
subdivergences are eliminated by the same counterterms as those which make finite the 
Bethe-Salpeter equation. This is most readily seen on the following form of the equation: 



9n(P) 



-\ I mQ)D\Q) + 



dm^ 



6{P-K)--r{P,K)D\K) 



(IV.18) 



dm'^ 2 Jq - - - - jk I 

obtained from Eq. pV.17j) by using Eqs. pi.2H) - pi.22j) . This expression of dU/dm^ in 
terms of F contains no vertex subdivergences. The remaining divergences can be absorbed 
in the mass counterterm, determined then by requiring that dU/dm"^ be finite. We shall 
illustrate this cancellation of subdivergences, by looking at a few iterations of Eq. ()IV.18|) . 
For the sake of simplicity, we use in all expressions the free propagator Dq, rather than the 
full propagator D (as we have already mentioned, this does not change the structure of the 
counterterms, only their numerical values; in particular the way the various subdivergences 
cancel is independent of this choice) . In zeroth order we have 



5ml 
m^ 



K,{P,Q)Dl{Q) 



(IV.19) 



where the symbol [■ ■ -j^.^ isolates the divergent part of the quantities inside the brakets. 
One can verify that Eq. pV.19|l leads to the expression pV.9|) of the mass counterterm 
obtained in perturbation theory. The first iteration gives 
5m\ 1 



m" 



4\P,Q)D'\Q) 



6ml 



(6{P-K)-ho{P,K)Dl{K)^,{W.20) 



which after using Eq. pV.12|l leads to 



Sml 



+ 



6ml 
m^ 



5X, + - I K,{P,K)Dt{K)ko{K,Q)]Dt,{Q) 



K 



jJ-UP,K)Dl{K). 



(IV.21) 



In contrast to the case of the 2-loop skeleton, where the multiple integrals factorize, here 
one cannot obtain an explicit expression. One can however easily recognize the role of 
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counterterms determined at a previous iteration (here Stuq), or those coming from the 
Bethe-Salpeter equation (here SXi), in the cancellation of subdivergences. Observe in 
particular that the product of 6mQ by the finite part of the integral in the second line, 
which is P-dependent, cancels with a corresponding term in the first line which arises by 
keeping K fixed and taking Q to infinity. 

C. The determination of 6Z and the asymptotic form of 11 

In the discussion of the two-loop approximation in Sect. II, we have been able to sepa- 
rate the calculation of the vertex and mass counterterms 6X and from the resolution 
of the gap equation for the propagator. This has been possible since the divergences ap- 
pearing in the Bethe-Salpeter equation for F, and also in Eq. pil.25|) for dU/dm'^, depend 
only upon the leading behaviour of the propagator at large momenta, and this is not mod- 
ified by the resummation performed by the one- loop gap equation pil.2j) . Therefore, in 
order to compute the counterterms 6X and 5m^, it was enough to consider these equations 
with the full propagator replaced by the free one, Dq{P) = (P^ -|- m^)~^ (where the mass 
m did not play any role except that of an infrared regulator). 

But in the present, three-loop approximation, such simplifications no longer hold. In- 
deed, as revealed by the analysis of the first orders of perturbation theory, the insertion 
of the self-energy does change the asymptotic behavior of the propagator, so that this 
behavior keeps changing as one proceeds through the iterations. Thus, if one were to 
compute the propagator and the various counterterms through iterations, both the gap 
equation and the equations for 6X and (5m^ would have to be solved simultaneously. Still, 
as we explain now, it is nevertheless possible to disentangle the calculation of the coun- 
terterms from that of the full solution to the gap equation. The key observation is that 
the counterterms are sensitive only to the leading asymptotic behavior of the propagator. 
This is clear in the previous perturbative analysis, and in the general case it results from 
the fact that Eqs. ()II.17|) and pV.18|) involve only logarithmic divergences. As for 6Z, it 
is determined by the global divergence of the self-energy diagrams obtained when all the 
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internal loop momenta are simultaneously sent to infinity (after all subdivergences have 
been removed): this divergence is proportional to , and is independent of the mass, or 
any other subleading contribution to the propagator at large momenta. 



Now, according to Weinberg's theorem 



15 



2J| , one expects that for very large momenta 



K ^ m, Il{K)/K^ grows, at most, as a (slowly varying) function of InK^. It is useful 
to introduce the special notation 112 (K) for the leading asymptotic contribution: by 
assumption (11 — 112)/-^'^ ^ as oo. To compute Il2{K), we need the limit of 

Eq. ()IV.2|) for very high (and zero temperature). Since the mass plays no role in this 
limit, one could set m = 0, and we shall indeed do that eventually. However in practice, 
it is sometimes useful to keep a non-zero mass to avoid infrared complications, so let 
us consider first the case m ^ 0. The corresponding gap equation involves all types of 
divergences, but, for the purpose of computing 112, it can here be renormalized by simple 
subtractions. That is, 112 (if) can be obtained from the large K behavior of the solution 
n to the following, auxiliary, gap equation: 

U{K) = S{K') - S{K' = -m') - + "^')^|^,_^,- (IV.22) 

where S[K) denotes the r.h.s. of Eq. pV.2|) . This equation has been used in Ref. [3| 
to compute the self-energy in the vacuum. It is finite indeed, as it can be checked by 
taking two derivatives with respect to within the integrand, and noticing that the 
result is finite. It should be emphasized that this procedure holds because, in the present 
approximation, the vertex counterterm (5A enters only a local self-energy diagram (the 
tadpole in the r.h.s. of Eq. pV.2j) ). and thus it can be simulated by a local, mass, 
subtraction. In any case, Eq. pV.22|) is correct only for the calculation of 112, and not for 
the subleading contribution IIq = 11 — 112. 

In what follows we shall compute 112 from the gap equation with m = 0, using dimen- 
sional regularization. Then, both the mass divergences and the vertex subdivergences 
automatically cancel out, and the remaining field divergences can be eliminated by the 
counterterm 6ZK^. The equation determining 112 reads then 

U2{K) = S{K) + K'^SZ (IV.23) 
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where this time S{K) denotes the r.h.s. of Eq. ()IV.2|) evaluated with propagator 
D2^(K) = + 112 (-fT). Clearly, the finite form of the equation, and thus also the 
ensuing solution for 112, will depend upon our choice of a subtraction scheme. For in- 
stance, one can choose 5Z = —{dS/dK'^)\K2^_^2, in which case the asymptotic behaviour 
of the propagator will depend upon the subtraction scale /x. Since, as we have seen be- 
fore, the counterterms 6X and (including their divergent parts) are sensitive to the 
asymptotic behavior of the propagator, it follows that these counterterms will end up 
being dependent upon the renormalization scheme for 6Z. 

We can now summarize the analysis of the last two subsections as follows: i) The 
values of the counterterms are determined by the dominant piece of the self-energy at 
large K. ii) They can therefore be determined from the solution of an auxiliary gap 
equation which coincides at asymptotically large K with the original gap equation, iii) 
The coupling renormalization involves two types of counterterms: those which renormalize 
the kernel A of the Bethe-Salpeter equation, and those which render finite the iterations of 
the Bethe-Salpeter equation. The counterterms of the first kind can be summed up as the 
counterterm 5\q. The counterterms of the second kind can be computed in terms of the 
solution to the Bethe-Salpeter equation, iv) The mass counterterms are determined by the 
equation for 811/ dm? where the vertex sub divergences are eliminated by the counterterms 
which render finite the Bethe-Salpeter equation. 

D. The gap equation: algebraic renormalization 

The verification that, with the counterterms determined as indicated above, the gap 
equation is indeed finite, goes through similar steps as in subsection 1111 PI The main 
complication here comes from the need to deal with the modification of the asymptotic 
part of the propagator resulting from the insertion of the momentum dependent self- 
energy. We shall handle this by considering an auxiliary gap equation, as indicated in the 
previous subsection. 
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We write the full propagator as 

where Il{K) is the renormalized self-energy. Then we set 

n = Ho + Ha (IV.25) 

where 112 is the solution of an auxiliary gap equation which coincides asymptotically with 
the original gap equation. As indicated in the previous subsection, 112 niay be obtained, 
for instance, as the solution of the gap equation with m = 0, that is: 

MK) = JpJq D2o{P)D2o{Q)D2o{K + P + Q) + K^5Z, (IV.26) 

where D2o{PY^ = + n2(P). Note that, thus defined, 02(0) = 0. We denote by Uo{K) 
the difference between the true solution, 11, and 112. It follows from our assumptions that: 

lim = 0, (IV.27) 

K^oo K 

SO that IIq{K) grows at most logarithmically at large K. There is of course some arbi- 
trariness in the decomposition pV.25jl . related to the specific choice one makes for 112. 
The important point is that the asymptotic behaviour of 112 (-ft') at large K coincides with 
that of Il{K). Indeed, from the previous discussions, we know that the counterterms 
depend only on 112 and are insensitive to the precise separation between 112 and IIo. 
The propagator 

D2{K) = {K^ + + Hs)-' (IV.28) 

will play the same role here as Dq{K) in the previous section, i.e., we shall expand with 
respect to it. We include im? within D2 in order to avoid possible infrared problems 
in this expansion (the massless case is discussed in Sec. II V F|) . We call r2 the solution 
of the Bethe-Salpeter equation where the intermediate propagators are D2 propagators 
and similarly the kernel A2 is evaluated with D2 propagators. We shall assume that 
the counterterms 5\ and 5m? are determined from this equation and the corresponding 
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one for dll/dm?. That is, the solution of the auxihary gap equation provides a finite 
112 and well defined values for the counterterms 5\ and 5m?. We verify now that 
these counterterms render the initial gap equation finite. The strategy will be to use 
the auxiliary gap equation to obtain D2 and then determine the counterterms. These 
counterterms will then be used in order to eliminate the divergences of the gap equation, 
leading to a finite equation for IIo. 

It is useful first to introduce a special notation for the self-energy calculated with 
the propagator D2- We set 112 (-ft') = n[D2], where n[-D] denotes the right hand side of 
Eq. (fTOjl : that is: 

ti2{K) = -^ [ [ D2{P)D2{Q)D2{K + P + Q)+^-^ [ D^iP) + Sm" + KHZ. (IV.29) 
Q JpJq 2 Jp 

Note that if the counterterms SX and Sm'^ in this equation are limited to their leading 
order values, 6\o and 6ml respectively (calculated with propagator D2), then 112 (-ftT) is 
finite. 

We now expand the propagator D around D2 

D = D2 + 6D 6D = D2 [-Ilo]D2 + Dr (IV.30) 

where we have singled out in SD the term linear in Ho. The usefulness of this expansion 
is that each extra power of Ho decreases the degree of divergence of the diagram in which 
the corresponding propagator is inserted. This property will allow us to isolate the vertex 
sub divergences in very much the same way as in the two-loop example. In analogy to 
what we did in Sect. IIII Dl we define 

flo{K) = ^J^A2{K,P)6D{P). (IV.31) 

Now we note that the gap equation can be written as 11 = Il[D] where D is expanded 
as in Eq. pV.30|l . and may thus be considered a function of Hq. Using the definitions 
()IV.29|) and pV.31|) above, we can put the gap equation in the form 

n = 62 + Ho + fir- (IV.32) 
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where 11 in the left hand side is to be understood as 112 + Hq (so that the equation is 
efffectively an equation for Ho). The quantity Ilr{K) is not only finite, but decreases as 
1 / at large K. This follows from Weinberg's theorem and the fact that 11^ contains at 
least two IIo insertions and is therefore of degree of divergence —2. 
Next we note that 

no(K) - no(0) = i ^ [A2{K, P) - A2(0, P)] 6D{P) (IV.33) 

is finite (cf. after Eq. (I1V.14|) ). Since both flo{K) -flo{0) and fl2{K) -fl2i0) are finite, so 
is therefore Il{K) —11(0). It follows that the ultraviolet divergences are entirely contained 
in n(0) = 112(0) + 60(0) + 11^(0) where both 112(0) and 60(0) are a priori divergent. 

For no(0) we proceed as in Sect. IIIIDl and express A2(0, P) in the defining equation for 
no(0) (Eq. (I1V.31I1 with K = 0) in terms of 12(0, P) (the solution to the Bethe-Salpeter 
equation with propagator D2). We get: 

no(0) = i f A2{0, P)SD{P) 



12(0, P) [6D{P) + Dl{P)flo{P) 
j^^2iO,P) [DriP) + Dl{P){IloiP) -UoiP))] (IV.34) 



i^|r2(0,P) + ^ J T2iO,Q)DliQ)A2iQ,P)^6DiP) 

1 
2 
1 
2 



For n solution of the gap equation, we can use Eq. pV.32j) to replace ITq by 11 — 112 — 
in the r.h.s. of this equation, and verify that the divergent terms linear in IIq (= 11 — 112) 
cancel, as anticipated. We are left with: 

no(o) = IJ^ 12(0, p) {DriP) - n.(P)z^2'(^)} 

+ ^ ^ r2(o, p) {n2(P) - n2iP)}DliP). (iv.35) 

The first line is finite, as it can be easily checked by using Dr{P) ~ 1/P^, n,.(P) ~ 1/P^, 
and the fact that r2(0,P) behaves like InP^. 
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The remaining divergences in 11(0) are therefore those of the following expression (we 
add 62(0) = 0): 

^ 12(0, P) {n2(P) -fl2(P) } DliP) + 62(0). (IV.36) 

We verify now that the counterterm Sm"^ absorbs the potential divergence. To do so, 
observe first that this expression vanishes in the limit m — * (then 112 (-P) = n2(P) and 
112(0) = 0). Let us then expand 112 — 112 = n[D2o] — n[Z}2] around D2 by using : 

where Ds{P) ~ 1/P^ at large P. We get 

U2{P)-fl2{P) = -Sm' + lm' [ A2{P,Q)Dl{Q) + U,{P), (IV.38) 

^ JQ 

where Ils{P) is finite and goes as at large P. Inserting this in Eq. ()1V.36|) . and 

keeping only the terms which are not manifestly finite, we obtain: 

T X { ^ - } ^'^^^ 

w|i-^r2(o,p)D2'(p)| 

(IV.39) 

One recognizes in this expression the right hand side of Eq. ()1V.18|) which specifies the 
value of the mass counterterm. It is therefore finite. 



E. Finite temperature 

It is straightforward to extend the diagrammatic analysis presented earlier in this 
section to the case of finite temperature. Doing so, one easily verifies that, iteration by 
iteration, all contributions involving the temperature are finite. Potential divergences 
coming from counterterms multiplying finite temperature contributions cancel out when 
sub divergences are eliminated. Furthermore, since the counterterms depend only on the 
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asymptotic part of the propagator, which is unaffected by finite temperature effects, they 
are manifestly independent of the temperature. (The fact that 112 (-P) is not modified 
by the temperature follows from Weinberg's theorem 24 1, which ensures that the 



asymptotic behaviour ~ P^F(lnP^) of the self-energy comes from vacuum diagrams in 
which the internal momenta in all the internal lines are large and proportional to P.) 

But although this diagrammatic analysis reveals that no difficulty of principle should 
be expected, in actual calculations the temperature dependent contributions may become 
mixed in a subtle way with vacuum, divergent contributions, and some work may be 
required to disentangle them. To be specific, we consider here the calculation in the 
imaginary time formalism. It is possible to organize the summation over the Matsubara 
frequencies in such a way as to separate the explicit thermal dependence in terms of 
integrals with 0, 1 and 2 thermal factors: 

U{K) = U^K) + U,^{K) + n,„(K). (IV.40) 

The first contribution, Ilo,-,{K) is exactly the same functional of D as at zero temperature, 
Eq. pV.2|) . Since D depends on the temperature (because 11 does), Ilon{K) depends 
implicitly on the temperature. The term with one statistical factor, Ilin{K), can be 
written as: 

U,,XK) = l [ cx{P)A{K,P) (IV.41) 



2jp 

where A is the same functional of D, Eq. pV.3|) . as at zero temperature and depends 
only implicitly on the temperature through D. We have set cr(P) = p{po,p)n\pg\ep^, where 
p{po,p) is the spectral function for the propagator D: 

D(.,p)=r^E^, (IV.42) 



n. 



p — l/(e^/^ — 1), and Sp^ = ±1 depending on the sign of Pq. Finally, we have introduced 



the notation P = {po,p), to be distinguished from P = {ipQ,p). The 4-point functions 
such as A2{K, P) are analytically continued in the complex plane from the complex value 
ipo where they are well defined (note that since A{K, P) is a function of K + P, it remains 
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well defined when P is continued to P). Finally the term with two statistical factors, 
Han (-ft'), can be written as: 

U^K) = yJpJ^ a{P)a{Q)D{K -P-Q). (IV.43) 

As in the previous subsection, we write 11 = 112 + 110, with 112 determined by the 
zero-temperature auxiliary gap equation pV.26|) . and IIq{K) is at most logarithmic and 
carries the whole dependence upon T. For Iio,-,{K) we set again: 

no„(K) = ti2{K) + no(ir) + tir{K), (IV.44) 

where 112, fto) and 11^ are given by the same Euclidean integrals as in the previous sub- 
section, but contain implicit temperature dependence through the propagator D, except 
for 112 which depends only on D2 and remains therefore temperature independent. It is 
also convenient to write 

n,„(i^) = ni(K) + n3(i^), (iv.45) 

where 

MK) = IJmP)MK,P) (1V.46) 

is at most logarithmic at large K, and can therefore be considered as a modification of Hq. 
As for II3, it decreases faster than a logarithm and it will be considered as a modification 
of 11^. The same remark applies to Il2n{K) that we shall denote by 114 (i^'). 

The final expression for the gap equation takes then a form similar to that at zero 
temperature, and may be written as 

n = n2 + (Ho + Hi) + (fir + 63 + 1I4) = 62 + n[, + n^, (iv.47) 

with IIq = Hq + Hi and 11^ = 11^ + 113 + 114. Applying the strategy of the previous 
subsection to eliminate the vertex sub divergences in nQ(0), one gets easily: 

^oW = lJ^DriP)T2iO,P) + \J^aiP)T,iO,P) 

~li r2(o, p)Dl{P) {u,{P) - ii,{P) - K{P)] . 

(1V.48) 
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This equation generalize Eq. ()IV.35|) . Its first line is finite. The divergences in the second 
line are temperature independent. When combined with those of 112(0), they become 
identical to those of the zero temperature case, and are absorbed in 6m?. 

F. Massless case 

The massless case is interesting only at non-zero temperature. (Indeed, if T = 0, 
then the gap equation with m = is simply the auxiliary gap equation that we have 
introduced before, cf. Eq. ()1V.26|) . and whose solution is 112.) Consider then the finite- 
temperature case, where 11 admits the decomposition shown in Eq. pV.40|) . When m = 0, 
the subsequent manipulations of the term with zero thermal factors, Hon, deserve a special 
treatement: Indeed, if one proceeds as before, i.e., by developing D around D2 as in 
Eq. pV.30|) . one meets infrared problems: since D2 is now a massless propagator (this is 
the same as D20 introduced before), the Bethe-Salpeter equation for the 4-point function 
r2 is infrared divergent, and so are also 11,. and Dr used in pV.35|) . One can however 
follow the same strategy as in Sec. IIII El and exploit the fact that the finite temperature 
effects generate a natural infrared cutoff. Since the ultraviolet divergences are ultimately 
sensitive only to the asymptotic behaviour of the propagator, they will not be affected by 
the replacement of D2 by D at appropiate places in the previous expansions. Specifically, 
we shall exploit the fact that 112 is the solution to the massless gap equation at T = 0, 
namely (cf. Eq. ()1V.26|) ) : 112 = non[i52]- Then, after rewriting the full propagator as in 



one can develop no,jD2] around non[-D] without encountering infrared divergences: 



Eq. (ITTLMII : 



D2iP) = D{P) + D{P)Il^{P)D{P) - D[{P) = D- 6D' , 



(IV.49) 




(IV.50) 



or, equivalently (compare to Eq. pV.44j) ): 




(IV.51) 
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where Uq^ = Ilon[D]- Note that, as compared to Eq. ()IV.31|) . the development above 
involves the kernel A built with the full propagator, and similarly for 6D'. It is then 
possible to use the complete Bethe-Salpeter equation (with propagator D) in order to 
replace A in terms of the infrared finite 4-point function F. 

The subsequent manipulations proceed as before, in Sect. II V El except for the replace- 
ments r2 — * r and 6D 6D'. Since 112 = ^2 when m = 0, this finally leads to the 
following, manifestly finite, expression for IIq (compare to Eq. pV.48|) ): 

no(^) = IJ^ r(K, P) [d^P] - Il',iP)D\P)] +^J^ a{P)T{K, P) . (1V.52) 

Since, however, the Bethe-Salpeter equation for F involves the full propagator D, and 
thus is sensitive to the complete self-energy 11 = 112 + IIq + 11^, it follows that, unlike in 
the massive case, the Bethe-Salpeter and the (full) gap equation are now coupled, and 
must be solved simultaneously. 



V. GENERALIZATION TO HIGHER-LOOP ORDER 

In order to generalize the previous considerations to higher-loop order skeletons, we 
note first that many of the results that have been obtained in the previous section are 
independent of the precise structure of the skeleton considered. In particular, the equa- 
tions which determine the vertex and mass counterterms, i.e. the Bethe-Salpeter equation 
flll.34|) and the equation ()1V.18|) for dU/dm^, have been used in the previous section in 
their general forms, without further specifying the diagrammatic structure of $, 11, or 
A. Similarly, the manipulations performed on these equations, like the expansion of the 
gap equation in Eq. ()1V.32|) . or the expression pV.31|) for IIq{K), are generic as well, 
and would hold for any choice of the skeletons in $[-D]. There are new features however 
that need to be discussed. They are already present in the 4-loop approximation to ^[D], 
and we shall illustrate our discussion with this example. The corresponding diagrams are 
displayed in Fig. |221 

As we have seen in the discussion of the 3-loop approximation in Sect. IV, skeleton 
diagrams with different number of loops get mixed under renormalization: the two-loop 
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skeleton is indeed needed as a counterterm to eliminate vertex subdivergences of the three- 
loop skeleton. We shall call primary skeletons the skeletons which are selected in the 
^-derivable approximation that one considers, and counterterm skeletons those skeletons 
which carry the counterterms needed to remove the vertex subdivergences of the primary 
skeletons. To be specific, let us order the skeleton diagams according to their number of 
loops, and consider the approximation in which the primary skeletons are (A^ + l)-loop 
diagrams (so that the corresponding diagrams for 11 involve N loops). The removal of 
divergences in the gap equation for 11 proceeds in three steps: 

i) In the first step, one determines the 4-point function A from 11 by using Eq. pi.l6|) . 
and one constructs the vertex counterterms which render A finite. Since obtained by 
opening one line in 11, the 4-point function A is given by primary diagrams with — 1 
loops. In general, these diagrams will involve vertex subdivergences. The corresponding 
counterterms are constructed as in ordinary perturbation theory, with however the prop- 
agator D associated to the lines of the diagrams; that is, they are the divergent parts 
of the loop integrals associated to the respective (sub)graphs evaluated with propagator 
D. These vertex counterterms are to be inserted in self-energy diagrams with L < N 
loops (the "countertem skeletons"), while the primary skeletons are calculated with the 
renormalized coupling at their vertices. 
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FIG. 23: The 2-loop contribution to A. 

In the previous examples, we have met one counterterm of this type, namely the coun- 
terterm 6Xq required to remove the global divergence of the one-loop contribution to A. 
For the 3-loop primary self-energy in Fig. |221b, the corresponding, 2-loop, contributions 
to A are displayed in Fig. 1221 and their various divergences are exhibited in Fig. |23 
These are both one-loop and 2-loop vertex counterterms, which are included in the ver- 
tices of the 2-loop and, respectively, one-loop self-energy graphs, as shown in Fig. [23 
Note that the one-loop counterterms, which cancel subdivergences of A, involve all the 
possible channels. By contrast, the 2-loop counterterm, which is a global divergence of 
this particular contribution to A, involves only the s and u channels (so like 6Xo in the 
example of Sect. IV). Clearly, this is related to the irreducibility of A in one particular 
channel, the t-channel. 




FIG. 24: The subgraphs enclosed in the dashed line boxes correspond to the divergences in A, 
and also to the vertex counterterms shown in Fig. 1251 The contribution to 5Xii^ represented by 
the diagram on the left corresponds to the s and u channels, while that on the right corresponds 
to the t channel. 
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X) 




(a) 



(b) 



FIG. 25: The counterterm diagrams necessary to remove vertex sub divergences from the 3-loop 
self-energy in Fig. 1221 b. 

This example can be generalized to all orders: For a A^-loop primary self-energy, the 
vertex counterterms associated with subdivergences in A include all the possible channels, 
and appear in the vertices of self-energy skeletons with 2 < L < — 1 loops (in such 
a way that L + I = N , where I is the loop number of the vertex counterterm). As for 
the global divergence of A, this gives a (A^ — l)-loop counterterm which excludes all the 
topologies which are 2-particle-reducible in the t channel. This counterterm is included 
in the one-loop self-energy diagram (the "tadpole"). 

Note that, as a result of this first step, the counterterms needed to renormalize A are 
not explicitly computed, but rather expressed as integrals which involve the (yet unknown) 
propagator D. However, the divergent parts of these integrals are sensitive only to the 
leading high momentum behaviour of the propagator, i.e., to Similarly, A (the sum 
of A plus the vertex counterterms constructed at this step) is a finite functional of D, or 



ii ) The second step refers to the renormalization and the calculation of the asymptot- 
ically equivalent propagator D2. To that aim, one proceeds as in the 3-loop example in 
Sec. IIV CI — that is, one considers the self-consistent gap equation with m = — except 
that, in the general case, one has to include in the r.h.s. of the equation for 11 all the 
counterterm diagrams determined in step i). Then, within dimensional regularisation, 
the functional n[D2] (the r.h.s. of the gap equation) has neither vertex, nor mass, diver- 
gences, but only a global field divergence, proportional to K^, that can be absorbed in 



D2. 
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the counterterm SZK'^. After this divergence is subtracted, the ensuing, finite, equation, 
can be solved for 112. This gives the propagator D2, which then can be used to exphcitly 
compute A2. 

Hi) In the third step, the (finite) kernel A2 together with the propagator D2 are used 
to construct the Bethe-Salpeter equation, together with the equation for 611/ Sm^. By 
renormalizing these equations, we obtain the countertermes 6X and (5m^ required to elim- 
inate the residual vertex sub divergences and, respectively, the mass divergences from the 
gap equation. Note that, in contrast to the vertex counterterms for A (cf. the first step), 
which involve, at most, A^ — 1 loops, those obtained by iterating the Bethe-Salpeter equa- 
tion will receive contributions with arbitrarily many loops (and similarly for 5m?). At this 
stage, one can use the algebraic derivation in Sect. IV. E to verify that the gap equation 
is a finite equation for Ho = 11 — 112. 

The arguments in Sect. IV. E rely on two properties which remain true in any order in 
the skeleton expansion, as we shall see. These properties are: 

A) A(P, K)—A{P, 0), considered as a function of P for fixed K, is of degree of divergence 
less than zero. 

B) The term 11^ (A') in the r.h.s. of Eq. pV.32|) is finite and decreases like at 
large K. 

The first property, on the behaviour of A(P, K) — A{P, 0) at large P for fixed K, follows 
from the 2-particle irreducibility of A, as we explain now. The different contributions to 
the logarithmic behaviour of the 4-point function at fixed K and large P are, according 
to Weinberg's theorem, associated to 4-point subgraphs attached to the two external lines 
carrying P. Recall that, in Weinberg's terminology, a subgraph is any set of lines choosen 
in the initial graph in such a way that there is no vertex attached to only one line of 
the subgraph jisi 2^ . In Fig. |2ni we give examples of 4-point subgraphs attached to 
P for a graph contributing to r{P,K). We also give an example of a set of lines which 
is not a subgraph. For both r{P,K) and A{P,K), there is always a trivial four point 
subgraph attached to P : the graph itself. Clearly the leading behaviour at large P 
for this contribution is independent of K, either for r{P,K) or A{P,K). The reason is 
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K K K K K K 




P P P P P P 



FIG. 26: Examples and counter-example of 4-point subgraphs attached to P. (a) and (b) are 
4-point subgraphs contributing to the dominant behaviour of T{P, K) at large P and fixed K. 
This leading contribution is independent of K in the case of (a) but the subgraph (b) may give 
a X-dependent contribution, (c) is not a subgraph. 

simply that for this subgraph, K appears in propagators where, by asumption, the total 
momenta goes to infinity, while K is kept finite. Of course there are other logarithmic 
contributions. The second diagram in Fig. |2ni gives an example of such a contribution in 
which the leading contribution at large P does depend on K. But it is easily verified that 
this particular 2-particle reducible topology is the only one allowing for a i^-dependence. 
Since A(P, K) is irreducible, its asymptotic behaviour at large P is independent of K. 
Property A) then follows. 

For the second property, we note that 11,. is obtained by differentiating twice a skeleton 
graph for the self-energy. Thus 11^ has the structure of a 6-point function which is closed 
into a two-point function by two 5D^ each of degree of divergence —4. The internal six- 
point subgraph itself has degree of divergence —2, and since all (vertex) sub divergences 
supposed to be eliminated, this is strictly finite. When the two effective propagators 5D 
are closed on it, the UV finiteness is preserved because one of the two following situations 
occurs: 1) The two propagators appear together in the same loop; that loop is certainly 
finite. 2) If a single effective propagator appear in some loop, that loop is certainly not a 
tadpole (since the whole graph has the topology of a skeleton), and thus contains at least 
one additional propagator belongning to the six-point vertex subgraph; that loop is finite 
too. Since the overall graph is finite, simple power counting together with Weinberg's 
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theorem ensures that, at large K, tlr{K) ~ 1/K'^. 

To conclude, let us briefly discuss the scale dependence of the renormahzed couphng 
constant which emerges from the renormalization program described above. (A similar 
discussion applies to the /i-dependence of 112 and of the renormalizcd mass parameter 
m.) Consider the A^-loop approximation to the self-consistent self-energy which includes 
all the primary skeletons with L < N loops, together with the corresponding counterterm 
skeletons. Given the similitude between step i) above and ordinary perturbation theory, 
it is clear that renormahzed vertices belonging to diagrams with different numbers of loops 
are dressed differently, and thus run differently with /x. The respective /i-dcpendences 
are non-perturbative, since contributions of all orders in A are included (because the 
corresponding counterterms depend upon 112, and also because of the iterations performed 
by the Bethe-Salpeter equation in the case of the coupling constant which enters the 
tadpole). However, when restricted to an expansion in powers of A, the //-dependences of 
the vertices are those expected from the loop order of the included diagrams. In particular, 
the couplings in the vertices of the A'"-loop primary diagrams are not dressed at all, so they 
are independent of /i, those which enter the {N — l)-loop primary diagrams are running 
according to the one-loop beta function, etc. The coupling constant in the one-loop self- 
energy is exceptional, in that it exhibits the most developped running (perturbatively 
correct up to loop order — 1). 

The one-loop example discussed in Sect. Ill (see, especially. Sect. III.D) suggests that 
it should be possible to exploit the non-perturbative character of the //-dependence of 
the vertices in order to minimize the scale dependence of the self-consistent solutions for 
r and n. In that one-loop case, we have in fact been able to verify explicitely that both 
r and n are strictly independent of /t. In the general case, the corresponding analysis 
is considerably more involved, and will not be pursued here. We shall only emphasize 
a specific feature of the ^-derivable approximations, which is that vertices in different 
diagrams run differently with /i. This should be taken into account when trying to reduce 
(and possibly eliminate) the scale dependence of the physical observables. 
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